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While the electromagnetic force is microscopically simply the Lorentz force, its macroscopic form 
is more complicated, and given by expressions such as the Maxwell stress tensor and the Kelvin 
force. Their derivation is fairly opaque, at times even confusing, and their range of validity all but 
a well kept secret. These circumstances unnecessarily reduce the usefulness and trustworthiness of 
some key quantities in macroscopic electrodynamics. 

This article presents a thorough yet pedagogical derivation of the Maxwell stress tensor and 
electromagnetic force in condensed media. It starts from universally accepted inputs: conservation 
laws, thermodynamics and the Maxwell equations. Simplifications are considered for various limits, 
especially the equilibrium, with a range of validity assigned to each expression. Some widespread 
misconceptions are scrutinized, and hidden ambiguities in popular notations revealed. 

A number of phenomena typical of strongly polarizable systems, especially ferrofluid, are then 
considered. In addition to enhancing the appreciation of these systems, it helps to solidify the 
grasp of the introduced concepts and derived formulas, and it demonstrates the ease with which the 
Maxwell stress tensor can be handled, inviting theorists and experimentalists alike to embrace this 
useful quantity. 
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I. INTRODUCTION 



A. The Coarse-Grained Force 
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(f L ) = (p e E) + (j e xB>. (1) 

Unfortunately, this is of little practical value, as we nei- 
ther have, nor indeed arc interested in, the detailed mi- 
croscopic information of the fields p e ,j e , E, B. On the 
other hand, although we do have the knowledge of the 
coarse-grained fields (p e ), (j e ), (E) and (B), or may ob- 
tain them via the macroscopic Maxwell equations, it is 
clear that the difference between the above true expres- 
sion and the "fake" Lorentz force, 

(Pe)(E) + (j e ) x (B), (2) 

may become quite important: Think of an iron nail in 
a magnetic field, which is obviously subject to an elec- 
tromagnetic force, although the "fake" force vanishes, as 
( Pe ),(j e )^0. 

The difference between the true and fake Lorentz force 
is frequently taken to be the Kelvin force, being PjV(-Bj) 
in the electric case [EJ , where Pj denotes the polarization, 
with summation over repeated indices implied. Orig- 
inally, the formula was microscopically derived: First 
write down the force exerted by an external field on a 
dipole, eg by taking the derivative of the electrostatic en- 
ergy with respect to the position of the dipole for given 
field. Then, assuming that the dipoles are too far apart 
to interact, and to feed back to the given field, take the 
total force density as simply the sum of all forces ex- 
erted on the dipoles in a unit volume. Clearly, this result 
neglects nonlinear terms in Pi (which account for inter- 
action) and is valid only in the dilute limit. This implies 
small polarization and therefore weak force, and one is 
free to add a second order term, writing the Kelvin force 
as 

P l V{{E l ) + Pi/eo) = PtVDt/eo. (3) 

There is a more sophisticated macroscopic derivation of 
the Kelvin force that we shall discuss next. Stating that 
the same constraint applies here, and one may still write 
the Kelvin force either as Pi'V(Ei) or as pV-Dj/eo is 
bound to raise a few eyebrows. But as will be shown 
in this review, it is really a rather straight-forward im- 
plication of the assumptions made during its derivation. 
By the same token, for linear constitutive relation, the 
Kelvin force is only valid for small susceptibility, ^€1, 
to linear order in \- 

The same holds for the magnetic case, where the corre- 
sponding force is much larger in praxis (as it is not sub- 
ject to ionization or residual conductivity). The Kelvin 
force expression is again valid only when weak, to linear 
order in the magnetization M <C (B), or the magnetic 
susceptibility x m ^ 1- And one may take its expression 
either in the usual form |5| , as 

M t V{{B t ) - p Mi) ee poMiVHi, (4) 



or equally justified, as MjV(-Bj). Consequently and re- 
markably, the Kelvin force is not usually a valid formula 
in ferrofluids, where frequently M s=s H and x m ~ 1 ~ 
though it is widely employed there. 

The nowadays generally accepted method to derive the 
electromagnetic, or "ponderomotive" force ff , cf §15 of 
Electrodynamics in Continuous Media by Landau and 
Lifshitz H , is a definit improvement over the microscopic 
one given above. It is completely macroscopic and starts 
from identifying ff as 

ff = -VMT -P(T,p)5 ij ], (5) 

where —11*°' denotes the macroscopic Maxwell stress ten- 
sor, and P{T, p) the zero-field pressure - the pressure that 
is there for given temperature and density, but without 
the external field. The ponderomotive force ff reduces 
to the Kelvin force, if the polarization and magnetiza- 
tion are (assuming the dilute limit) proportional to the 
density p. 

The idea behind the identification seems attractively 
simple: Since — VjLT* * is the total force on a volume ele- 
ment, and — ViP the force at zero field, their difference 
ff must be the force due to the external field. Unfortu- 
nately, as we shall show in this review, the resultant force 
is rather ambiguous. The reason lies with the zero-field 
pressure P(T, p), which is not a unique quantity, because 
its value depends on the chosen variables, eg P(T, p) ^ 
P{s,p), where s is the entropy density. The convention- 
ally defined Kelvin force is valid only under isothermal 
conditions with the density kept constant. Under adi- 
abatic conditions, with s given, the temperature T is 
known to change with the field, so will P(T,p), the al- 
leged zero-field pressure. Any field dependence implies 
that P(T, p) contributes to the electromagnetic force. 
Employing the force of Eq (H) alone then leads to wrong 
predictions. 

In physics, we tend to start from a robust base of broad 
applicability when studying a given situation, narrowing 
our focus along the way only when necessary. Clearly, 
while the universally valid Lorentz force is one such base, 
the fragile ponderomotive force, and even more specific, 
the Kelvin force are not. This makes it desirable to look 
for better alternatives, one of which is employing the 
Maxwell stress tensor directly. 

B. The Maxwell Stress Tensor 

Given in terms of thermodynamic variables, the Max- 
well tensor is an unequivocally macroscopic, coarse- 
grained quantity. It holds for all conceivable systems and 
is subject only to the validity of local equilibrium - a con- 
straint in frequency but not in the type and strength of 
interaction. Besides, more frequently than not, in under- 
standing the many facets of strongly polarizable systems, 
it is more convenient to employ the Maxwell tensor di- 
rectly, rather than to follow the detour via pressure and 
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the ponderomotive force ff — in spite of having to deal 
with quantities such as the chemical potential or the en- 
tropy. 

Given the Maxwell tensor's pivotal role for a solid un- 
derstanding of electromagnetism in condensed systems, 
(not least as the starting point for forces), it is espe- 
cially unfortunate that one is hard pressed to find a fully 
convincing derivation of its form. The classic reference 
is again to the otherwise excellent book by Landau and 
Lifshitz ||, the relevant §15 of which, however, leaves 
many readers unconvinced, even confused. The reasons 
for this are presented in the four sections below. 



Nevertheless, it is customary to take for granted that 
the form of the Maxwell tensor and the electromagnetic 
force thus derived also hold in dynamic situations, at fi- 
nite frequencies and nonvanishing liquid velocity. This 
may well be approximately true for some circumstances, 
but which and why? Should we not rather derive dynam- 
ically valid expressions for the Maxwell stress tensor and 
the electromagnetic force, and justify the employed ap- 
proximations by estimating the magnitude of the terms 
neglected? 



3. A Technical Flaw 



1. Three Obvious Objections 

The three more obvious reasons for the frustration of 
the readers are (i) the gradient terms, (ii) the magnetic 
field dependence, and (iii) the concentration: First, the 
Maxwell tensor is derived considering a capacitor. Its 
geometry is so simple that terms proportional to the gra- 
dients of temperature and field may have been missed. 
Although these are asserted to be negligible in a foot- 
note, one feels the need for convincing arguments. 

Second, the magnetic terms are obtained, in §35 of ||, 
by the replacement E — * H, D — > B. This would be jus- 
tified if the system under consideration does not contain 
sources, and the respective Maxwell equations are the 
same, VxE = 0, V D = versus VxH = 0, VB = 0. 
In the geometry used to derive the electric part of the 
stress tensor, however, the charge on the capacitor plates 
is essential for producing the electric field. Therefore, an 
independent consideration for the magnetic part is nec- 
essary, with a current carrying coil replacing the charge 
laden capacitor. In §4 of his widely read book on fer- 
rofluids H, Rosensweig points to this gap, and aims to 
fill it. Unfortunately, the calculation starts from a faulty 
assumption, discrediting the derivation, cf section HIE 2 . 

Finally, neither Landau and Lifshitz nor Rosensweig 
study the concentration dependence of the stress tensor, 
although it is an important thermodynamic and hydro- 
dynamic variable in ferrofluids, liable to undergo much 
greater variations than the density, and hence of experi- 
mental importance. 



2. The Range of Validity 

Next, the lack of any discussion on the range of va- 
lidity for the obtained expression further undermines the 
trust of the readers. The derivation is thermodynamic 
in nature and takes place in a stationary capacitor, so 
the result is apparently valid only in equilibrium and for 
vanishing liquid velocity. Consistent to this, the static 
Maxwell equations for the rest frame are repeatedly em- 
ployed when deriving the electromagnetic force from the 
Maxwell tensor. 



Then there is a technical flaw that actually renders the 
derivation mathematically invalid: There is no question 
that one can always obtain a tensorial quantity 11'°* by 
evaluating the scalar ED° ^rij while varying the unit vec- 
tors £i and rij. Taking each of the two unit vectors to 
successively point in all three directions, we obtain nine 
different results for 11*°* ^nj, each a component of 11*°*. 
In ||, rij is the surface normal, £j the direction of the 
virtual displacement, and the Maxwell tensor — II*°* is 
denoted as o~ij. 

With the electric field being the only preferred direc- 
tion, £j and Uj need to be both along and perpendicu- 
lar to the field. In §15 of ||, only the displacement of 
the capacitor plate is considered, implying that the sur- 
face normal rij is kept parallel to the E-field throughout. 
Therefore, if the direction of the electric field is taken as 
e x , only the components T1 XX7 Tl yx and Tl zx could have 
been evaluated. A footnote states: It is not important 
that in this derivation E is parallel to n ; since Uy can 
obviously depend only on the direction of E, not on that 
of n. This is incorrect, because although Ily does not 
depend on n, the scalar Hij^rij does: Take for instance 
the term DiEj, which is part of the Maxwell stress, and 
which indeed does not depend on n. Yet with D, E || e x , 
we have DiEjUj ^ for n || e x , and D. L Ejnj = if 
n _L e x . 



4- Conceptual Problems 

Finally, there are conceptual problems of more funda- 
mental nature. The Maxwell stress tensor is introduced 
in || by the terse sentences: It is well known that the 
forces acting on any finite volume in a body can be reduced 
to forces applied to the surface of that volume. This is 
a consequence of the law of conservation of momentum. 
The surface force density is subsequently taken as the 
stress tensor. 

Momentum conservation is without further elaboration 
a confusing argument, because the very idea of force is 
something that alters the momentum of particles. Being 
accelerated by external fields, the momentum of a lump of 
material is not conserved, and there is no reason why the 
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momentum density pv should satisfy a continuity equa- 
tion. On the other hand, what the electromagnetic field 
does is to exchange momentum with the material, so the 
total momentum of both remains conserved, which is, as 
we shall see, g tot — pv + ~ExH. This momentum density 
does satisfy a continuity equation, gl' 



= 0. 



(Even in the presence of fields that accelerate particles, 
the total system - including particles and the field pro- 
ducing charge carriers - is still translational invariant, 
with the associated total momentum a locally conserved 
quantity, see §75 of g, Ch 14 of §, and §.) The neg- 
ative of the associated momentum flux, is the 
properly defined, dynamically valid Maxwell tensor, ft 
holds its ground for finite frequencies, nonvanishing ve- 
locities of the medium, and including dissipations. Un- 
derstanding the conservation of gr tot = pv +E x H renders 
the derivation of the Maxwell tensor cogent and coher- 
ent, even though the second term is, under usual circum- 
stances in condensed matter, much smaller than the first. 



C. Content of this Article 

Our aim in this paper is a careful and rigorous deriva- 
tion of the Maxwell tensor and the dynamically valid elec- 
tromagnetic force, in a presentation that is detailed and 
pedagogical. We also consider applications, especially for 
experimentally relevant situations in ferrofluids, to give 
hands-on illustrations on how to employ these quantities. 

The paper may be divided into two parts, The first con- 
sists of Chapter || and concentrates on the fundamental 
concepts. The hydrodynamic equations for a dense sys- 
tem exposed to electromagnetic fields — including espe- 
cially the form of the Maxwell tensor — are derived here. 
The second part (Chapter 111 and IV) is more practically 
concerned. It starts from the Maxwell tensor, derives 
the electromagnetic force, considers different simplifica- 
tions, and studies a number of experimentally relevant 
situations. Although we believe that Chapter y provides 
useful and necessary insights, the second part could if 
necessary be read alone. 

Chapter [n| starts by differentiating two types of the- 
ories: First, the high- resolution theory for a low den- 
sity system, with at most one particle per infinitesimal 
volume element, or per grain (as in photographs); sec- 
ond, the low resolution theory for a high density sys- 
tem, with many particles per grain. The microscopic 
Maxwell equations and the Newtonian equation of mo- 
tion (including the Lorentz force) belong to the first type. 
There are no hidden charges, polarization or magnetiza- 
tion here, and we know the whereabouts of every single 
particle. The second type is represented by the macro- 
scopic Maxwell equations, or any thermodynamic and 
hydrodynamic theories. The question about the coarse- 
grained form of the electromagnetic force arises here. 

We shall consider three systems with increasing densi- 
ties, to be accounted for by theories of decreasing resolu- 
tion and increasing complexities: The first system (sec- 



tion [I A) is a dilute gas of charged particles, which is 
well accounted for by the Newtonian equation of motion 
and the micr osco pic Maxwell equations. The next sys- 
tem (section [IB ) is a slightly dissociated liquid of par- 
ticles possessing negligible electric and magnetic dipole 
moments. So there are many particles, but at most one 
charge carrier, per grain. The appropriate theory here 
is of a mixed type, being a combination of the hydrody- 
namic theory and the microscopic Maxwell theory. The 
former accounts for all particles, the latter accounts for 
the spatially slowly varying field, either externally im- 
posed or p rodu ced by the few charge carriers. Finally, 
in section II C, we consider a dense system, including 



dipole moments and hidden sources. It is to be accounted 
for by the genuinely low-resolution, hydrodynamic-type 
Maxwell theory that is our goal to derive and consider in 
details. 

In all these systems, there are a conserved energy u tot 
and a conserved momentum density, <?,- ot ^ pvi, both con- 
sist of material and field contributions. These we shall 
choose as two of our dynamic variables, with the equa- 
tions of motion u tot + VkQf* = and g\ ot + V fe n^* = 0, 
where —11*^ is the Maxwell stress tensor. We shall ob- 
tain explicit expressions for g\ ot and IB?* in all three sys- 
tems — whereby finding the expressions in the genuinely 
dense, third medium is rendered very much simpler, and 
more transparent, by the considerations in the previous 
two dilute systems. 



In Chapter III, starting from the Maxwell tensor, we 
shall first obtain the general expression for the electro- 
magnetic force density, valid for arbitrarily strong fields, 
moving medium and finite frequencies, given essentially 
by — (sVT + pV£), where £ is the chemical potential. 
Although this expression is already valid without field, 
it is the proper formula for the quantity usually written 
as the sum of the zero-field pressure gradient and the 
ponderomotive force. In contrast to the latter, however, 
— (sVT + pV£) is valid for any sets of variables, and for 
arbitrarily large polarization and magnetization. 

In equilibrium, since both the temperature T and 
the chemical potential £ are uniform, this force density 
is identically zero. Nevertheless, electric and magnetic 
fields remain operative, may lift a polarizable system off 
the ground, deform its shape and generally alter its hy- 
drostatics — all signs that forces of electromagnetic ori- 
gin are at work. In section [IIP , the responsible equi- 
librium force is identified as a surface force. Its form 
follows directly from the simplifications that occur for 
the M axwell stress if the system is in equilibrium (section 
[II C| ) . As most of the experiments considered in Chapter 
[V reflect varying appearances of the equilibrium force, 



this is a useful and concise result. 

A thermodynamic rederivation of the Maxwell stress 
tensor is given in section [II E , an approach which is sim- 
ilar in spirit to that of Landau and Lifshitz P, ^| , yet 
carefully freed of the flaws listed in section |IB| above. 
In comparison to the derivation given in Chapter 0, this 
here is narrower in focus and validity, but simpler. Chap- 
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ter ID ends with section [II F, which points out the in- 
herent ambiguity of the zero-field pressure P(T,p), as 
it is independent of external fields only if T and p are 
kept constant. We stress this point by repeating the 
separation of the zero-field pressure from the Maxwell 
stress for various combinations of variables, obtaining 
many forms of the ponderomotive force, all rather dif- 
ferent from each other. Furthermore, a warning is issued 
on the unsuspected, narrow range of validity for various 
expressions in the notation containing the zero-field pres- 
sure, especially the Maxwell stress in the approximation 
(P + ^H 2 )Sij — HiBj, and the Kelvin force. Neither is 
valid for large magnetization and magnetic susceptibility, 
M fa H and X m ~ 1- 

The very useful magnetic Bernoulli equation by Ro- 
sensweig Q is unfortunately also given in this notation. 
An formulation free of ambiguity is found, along with the 
generalization to include the variation of concentration. 



vertical, current-carrying wire that goes through a dish 
containing ferrofluid will drag up the ferrofluid along the 
wire. Our calculation in section |IV B includes the surface 



In Chapter IV , we study a number of phenomena typ- 



ical of strongly polarizable systems. Most of these ex- 
periments have already been given detailed theoretical 
analyses, especially by Rosensweig in his comprehensive 
book on ferrofluids g. The reason to repeat them here 
is threefold: First, to demonstrate that it is not nec- 
essary to introduce the zero field pressure, that experi- 
mental phenomena may be well and easily accounted for 
employing the Maxwell stress directly, using a generally 
valid, unambiguous notation. Second, we believe that our 
derivations are well streamlined and easy to follow. This 
is achieved by stressing the simplification of equilibrium, 
and by consistently using only the Maxwell stress and the 
associated boundary conditions. In contrast, confusingly 
many pressures are deemed necessary in H- Finally, we 
aim to amend some of the physics and especially include 
the concentration of magnetic particles as an indepen- 
dent variable. Being inhomogeneous where the field is, 
the concentration alters experimental outcomes and is 
something theoretical considerations need to include. 

Since mesoscopic particles diffuse slowly, the equilib- 
rium between concentration and field gradient needs time 
(up to days, even weeks) to establish. So for times much 
briefer, the concentration will remain constant and may 
indeed be taken as dependent. On the other hand, strong 
gradients, small geometry, large particles will all consid- 
erably shorten the equilibrating time (down to fractions 
of a minute) , then the spatially varying concentration will 
be an important feature in any system. Besides, long con- 
tacts with field gradients are also of interest (eg in mag- 
netic O-rings). Although the non-uniform equilibrium 
distribution of magnetic particles is not a well studied 
subject, their diffusion under enforced non-equilibrium 
situation is, see H, ||. 

Now a rundown of the considered experiments. In sec- 
tion IV A , we study the variation of density and concen- 
tration in the presence of inhomogeneous fields: While 
the variation of density is a measured phenomenon, that 
of the concentration is not - although its variation is 
much more pronounced, by many orders of magnitude. A 



tension and t he variation of concentration. The following 
section IV C considers how far an external field will raise 
the level of a fluid column, and what a pressure gauge 
actually measures within the f errofluid that is subject to 
external fields. Section IV D deals with rather similar 



physics and co nsider s magnetic O-rings and scrap sepa- 
ration. Section IV E contains a simple calculation on the 
elliptical deformations a droplet of ferrofluid undergoes 
along the external magnetic field. 



D. Two Caveats 

There are two points we shall not address in this arti- 
cle. First, we shall not consider any dissipative terms, 
although the framework presented here is well suited 
for including them. This was a difficult decision, as 
the dissipative terms account for irreversibility, without 
which no macroscopic, coarse-grained theory is complete, 
and because they give rise to additional electromagnetic 
forces of comparable magnitude and experimental rele- 
vance. Yet, including these terms would have consid- 
erably lengthened an already lengthy article, and fur- 
ther complicated its formulas and arguments. We pledge 
amendments for the future, and for now direct the impa- 
tient readers to the literature || |], ^| . 

Second, given the long and tortuous history of the 
struggle to come to terms with the macroscopic electro- 
magnetic force, time and again forcing us to back up from 
blind alleys, any attempt by us on a comprehensive cita- 
tion would bear historic rather than scientific interests. 
As we see it, the consistent thermodynamic treatment 
given first by Landau and Lifshitz ||, ancl later general- 
ized to include nonlinearity by Rosensweig and others Q , 
cf the excellent summary by Byrne JTo| | , is so much more 
superior than the many treatises preceding them, that 
this represents the only worthwhile starting point for any 
further considerations — and consequently these two are 
also the only ones subject to scrutiny and criticisms here. 



E. Notes on Units 

This paper will be in the SI units throughout, though 
with a little twist to render the display and manipulation 
of the formulas simple. We define, and from now on shall 
extensively use, the fields and sources 



H = Hy/jIE, B = B/ y/Jj^, g e = g e / ^JT , 
E = EV^, D = D/V^, q = q/^T , (6) 
P = D — E = P/^r , M = B H = MVT^, 

where (again starting from here) the quantities with hats 
are the usual ones, denominated in MKSA, or the SI- 
units. The advantage is, all new fields have the same 
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dimension, yJ/ ni 3 , and sensibly, H = B and D = E in 
vacuum. (Electric charge q has the dimension \J Jm, and 
its density [g e ] = y/j/ m 5 . ) Written in the new fields, 
all formulas are rid of the ubiquitous (and content-free) 
factors eo and To return to the more familiar MKSA- 
fields at any instance, simply employ Eqs (^). (These 
fields are sometimes referred to as being in the Heaviside- 
Lorentz, or the rational, units. But saying so immedi- 
ately leads to heated and pointless discussions. We stick 
to the SI units throughout and use the above fields only 
as a convenient shorthand that may be abandoned when- 
ever one chooses.) 

Some literature, notably Q, employs cgs, or Gauss- 
waits. To ease comparison with them, simply note that 
the Gauss-fields, with tilde, may be obtained any time 
via 



H = H/V4^, B = B/V47r, Q e 
E = E/V4^, D = D/V47r, j e 
P = V4^P, M 



47rj e , 
IttM, 



(7) 



II. HYDRODYNAMIC THEORIES 

A. Dilute Plasma 

We first consider a rarefied system of charged particles, 
and choose a resolution that is high enough that each 
grain (or infinitesimal volume element) contains one or 
less particles. (This theory is meant as a starting point, 
to clarify a few concepts important for the more complex 
theories of the following chapters. So we shall simply 
discard the possibility that even in a rarefied gas, two 
particles will occasionally come close to each other.) The 
microscopic Maxwell equations account for the time evo- 
lution of the electromagnetic field for given sources, 



V e 

cVxb 



; Pe, 
PeV, 



V-b = 0, 

b = — cV x e. 



(8) 
(9) 



while the feedback of the field on the motion of the 
sources is given by the Newtonian equation of motion, 



mv a = q(e + v a x b/c), 



(10) 



one for each particle a. (To emphasize the fact that we 
are here dealing with the high-resolution, microscopically 
accurate fields, these are denoted by the lower case letters 
e = e.y/eZ, b = b/ \fpo, while the coarse-grained fields 
below are denoted by the usual D, E, H, B.) Eqs (|,|0) 
represent a conceptually simple and complete theory, 
but it contains a notational inconsistency: The Maxwell 
equations are an Euler type theory, accounting for the 
time evolution of fields at a given point in space, while 
the Newtonian equation is of the Lagrange type, which 
concentrates on a given particle. (So the term p e v in 
Eq (J9J) denotes the electric current at a space point, while 
v a in Eq (|l^) is the velocity of a particle.) As only the 



Euler formulation lends itself to a scaling- up of the grains 
and a reduction of the resolution, we shall first find the 
Euler formulation for the Newtonian equation, before us- 
ing it to draw a number of conclusions useful for dense 
systems. 



1. Euler Version 0} Newtonian equation 

Since there is at most one particle per grain, of vol- 
ume Vq, we may identify the velocity, mass and charge 
of a volume element with that of the particle occupying 
it at a given instance, and take all three to be zero if 
there is no particle. As a result, we obtain three (highly 
discontinuous) fields: v a — > v(r, t), m/Vc — > p(r,i), 
q/Vq — ► p e (r,t). The many Newtonian equations of mo- 
tion then reduce to one field equation, 



pj^v = p[v + (v ■ V)v] = p e (e + v x b/c). 



(11) 



(The quantity Vt> is to be taken from the velocity of 
the same particle at two successive moments.) It is now 
essential to explicitly include the continuity equation, 



p+V-(pv) = 0, 



(12) 



that in the Lagrange version is implicit, nearly inciden- 
tally contained in the fact that one does not loose any 
of the many equations ([l0]) . [The continuity equation for 
the charge density p e is implied by the Maxwell equations 
Eq§!).] 

The material contributions to the energy and momen- 
tum density are, respectively, 



p(c 2 +v 2 /2), 



,M 



pv. 



(13) 



where the first comprises of the rest energy and the (non- 
relativistic) kinetic energy. Employing Eqs ( |H] , |l2"| ), we 
find 



•,M 



V Q 



M 



Pe(V 



■ M 



V k Il%=p e {e + vx b/c), 



q m = u M v, nr k = gr v k . 



2. The Field Contributions 



(14) 
(15) 
(16) 



The field contribution to the energy and momentum 
density are 



/ = i(e 2 



b 2 ), g F = e x b/c. 



From the Maxwell equations (||||) we deduce 

-p e (v ■ e), Q F = ce x b 

9. 

= (e 2 + b 2 - u F )5. 



Vfcnf fe = -p e (e + v x b/c)i, 



IT 



b,b 



iv k . 



(17) 

(18) 
(19) 
(20) 



Note the relationship g F — Q F /c 2 . This is far from 
accidental and derives from the symmetry of the rela- 
tivistic energy- momentum 4-tensor, H F a = Iljg a , because 
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F 

c 9k 



TT F 



f/c = n F fe . Less formally, g F = Q F /c 2 
may also be seen as the conservation of the field an- 
gular momentum for p e — > 0. The angular momentum 
density £ F = r x g F is a locally conserved quantity in 
neutral systems. Rewriting Eq (n9h (with p e — > 0) as 
<%f /dt + V m (e y fcr : ,n F m ) = e,jfen|~7we observe that the 
angular momentum if satisfies a continuity equation only 
if the stress tensor is symmetric. Although this argument 
seems to require merely the symmetry of the momen- 
tum 3-tensor, Ilf k — II F j, we know that a nonvanishing 
n F 4 — n F fc in one inertial system will foul up the symme- 
try of the 3-tensor in other systems, as the antisymmetric 
parts of any 4-tensors mix in a Lorentz transformation. 
Yet angular momentum is conserved in every inertial sys- 
tem. 

Since this reasoning is so general, it also holds for the 
material part, g M = Q M /c 2 . Hence the expression for the 
momentum density is actually g M = pv[l + v 2 /(2c 2 )], cf 
Eqs ([T^ , [l6| ), though we are quite justified to neglect the 
second term in the nonrelativistic limit. Later, when we 
have no prior knowledge of the form of the momentum 
density, we shall deduce it from that of the energy flux, 
as angular momentum is also conserved in dense systems. 

We register the fact that while the expression for the 
energy density u F = ^(e 2 + b 2 ) is a genuine input - in- 
dependent of, and in addition to, the Maxwell equations, 
the formula g F — ex b/c is not, since Q F is given once 

F ■ 

U IS. 



3. Energy and Momentum Conservation 

The preceding two sections allow the simple and note- 
worthy conclusion that our starting equations imply local 
conservation of total energy, momentum and angular mo- 
mentum even in the presence of charges, p e ^ 0. Taking 



u tot = U F + U M and g tot = g F + g M ; we find 



V-Q tot = 0, <7* ot + V fc n*f=0, (21) 



nS t = ng t = n& + nS, Q tot = q f + q m , (22) 



Q tot /c 2 



pv + exb/c. (23) 



These results have been collected here because local con- 
servation of these quantities is always true, independent 
of the above derivation tailored to a dilute and finely 
resolved system. So we may start from them as input 
next. 



B. Weakly Dissociated Gas 

Now we consider a dense macroscopic system in its hy- 
drodynamic regime: To the above gas of dilute charge 
carriers we add a dense system of neutral particles with 
vanishing electric and magnetic dipole moments. One 
example of this composite system is a slightly ionized 
gas with negligible electric and magnetic susceptibility. 
This is still a comparatively simple system, as the highly 



resolving, vacuum Maxwell equations (^jj^) remain valid 
- and with them every single of the equations of sec- 
tion II A 2, about the field contributio ns to energy and 
momentum. The equations in section II A 1, about the 



respective material contributions, must be modified, as 
these will now be accounted for by three hydrodynamic 
variables: the entropy density s, and the coarse-grained 
mass density (p) and velocity (v) - all smooth and slowly 
varying fields. (The coarse-graining brackets are dropped 
below to retain the simplicity of display.) 



1. The Material Contributions 

The hydrodynamic theory of the material part of our 
system is given by two continuity and two balance equa- 
tions, 



p + V ■ (pv) = 0, s + V • (sv) = 0, 
u M + V-Q M =p e (v-e), 
fff + V fe n^ = p e (e + «xb/c) i , 
QM = (u M + P)v, Tlf k = gfv k +P5 ik . 



(24) 
(25) 
(26) 
(27) 



For p e — > 0, these are the well known hydrodynamic 
equations of a neutral, iso trop ic liquid p"l| . For finite 
p e , the right sides of Eqs (^5[ ^) have the given form 
because the field source terms of Eqs (JT^, |l^) remain un- 
changed, and because summing up the respective right 
sides must yield nil, such that total energy and momen- 
tum are conserved, as given in Eqs@. (Since v actu- 
ally denotes (v) now, the above formulas presume the 
validity of e = cV x b — p e (v), rather than the origi- 
nal Eqs @. This is justified, however, because charge 
carriers do move with (v) in the absence of dissipation.) 



Comparing Q M and Tlfj here to those of Eqs ( |16[ ) , the 



only apparent difference is the appearance of the pressure 
P, clearly an expression of the present system being dense 
and interacting. However, there is more than meets the 
eyes, and the most important one is our ignorance about 
the explicit expression for the energy u M . (Also, being 
a hydrodynamic, coarse-grained theory, the equations of 
motion now contain dissipative terms to account for ir- 
reversibility and entropy production. As mentioned in 
the introduction, however, we shall neglect these in this 
article.) 

In the co-moving, local rest frame of the liquid, (de- 
noted by the subscript o,) the material energy density in 
local equilibrium is given by the thermodynamic expres- 



du™ = Tds + Cod/9, 



(28) 



which states the simple fact that is a function of 
s and p, with the temperature T and the chemical 
potential £o defined respectively as T = du^/ds and 
Co = du^f I dp. (The dimension of the chemical potential 
Co is energy /mass.) The explicit form of u^ 1 is usually 
unknown, it depends on microscopic specifics, and is in 
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general rather complicated. Yet, as we know from ther- 
modynamics, and as we shall see below, Eq ( p8| ) is as 
such already very useful. 

The second of Eqs ([l3|), g M — pv, remains an excellent 
approximation, because the energy flux Q M as given by 
Eq (E7|) is still dominated by the term pc 2 v (contained 
in vr^v), from the transport of rest energy. All the other 
terms are (in comparison) relativistically small and can 
be neglected. 

In nonrelativistic physics, it is customary to subtract 



d(c p) = c dp from both sides of Eq (|28j), eliminating the 
rest energy from dirff on the left, redefining the chemi- 
cal potential as £o — c 2 011 the right, and at the same 
time subtracting c 2 pv from the energy flux Q M , Eq (27). 
This represents a shift to a different but equivalent set 
of independently conserved quantities: from u M and p to 
(u M — pc 2 ) and p. However, this operation does alter the 
link between the momentum density and the energy flux, 
rendering it g M = Q M /c 2 + pv « pv. We shall follow 
this convention, so and £cb starting from this section, 
no longer include the rest energy and c 2 , respectively. 

To account for a liquid with a varying local veloc- 
ity v, we need to employ a globally valid "lab"-frame. 

pv 2 /2, 
So the 



,M 



The material energy is then given as u M — u, 
where d(\pv 2 ) = d(g M /2p) = v ■ dg M - \v 2 dp 
thermodynamic expression of Eq (Eq) becomes du M = 
Tds + £dp + v- dg M , with £ = £ —~^v 2 the chemical po- 
tential of the lab frame. Clearly, this expression signals 



M 



that u — u (s, p a ,g ) is also a function of g 

If we are dealing with a solution, a suspension, or any 
system with more than one conserved densities, we need 
to replace the term £dp with £ a dp a , introducing all con- 
served densities p a as independent thermodynamic vari- 
ables. 



du M = Tds + £ a dp a +v-dg 



M 



(29) 



where summation over a = 1, 2, ■ • • is implied. Being 
conserved, all p a obey continuity equations, so the first of 
Eqs (HJ) is to be replaced by p a + V- (p a v) — 0. Usually, 
it is more convenient to retain the total density p as a 
variable. For a binary mixture with the solute density 
pi and the solvent density p2, we therefore take p\ and 
p = p\ + pi as the variables, implying that £ Q dp Q in 
Eq (|29| ) is to be read as £dp +£idpi. 

The pressure, defined as the energy change if the vol- 
ume changes at constant entropy and mass, is related 
to the variables and conjugate variables of u M via the 
Duhem-Gibbs (or Euler) relation, 



P = -d(Jufd 3 r)/dV = -< + Ts + {i ) aPa 



,M 



Ts + £ a p a + v ■ g 



M 



(30) 



Combining Eq (E9) with (BOh yields 



VP = sVT + p a V£ a + gfVv 3 . (31) 



Inserting this and the first of Eq (M) in (f26|), we find 
Pi± v = Pe{e + v x b/c) 
-(sVT + Pa V£ Q + gfVv,), (32) 



an equation that will prove useful later. Note an inter- 
esting feature of the hydrodynamic equations of neutral 
systems, Eqs ( |2^p5| , |32] ) with p e — > 0: They are given in 
terms only of the quantities appearing in Eq (|29j): the 
energy w M , the thermodynamic variables s,p a ,g^, and 
the conjugate variables T,£ ai v. Without an explicit ex- 
pression for u M , the hydrodynamic theory is clearly nec- 
essarily written in these general and abstract quantities, 
and we may take this observation as an indication that 
the hydrodynamic theory contains no more input than 
conservation laws and thermodynamics, which is the ba- 
sic reason of its general validity. We shall return to, and 
build upon, this point in the next section. 

For a charged system, p e ^ 0, we need to include the 
Maxwell equations for a complete description. The 
independent variables now include e and b, each with 
an equation of motion. All equations also depend on p ei 
which however only stands for V ■ e and is not indepen- 
dent. 



2. "2-field" Theory 

For our purpose of deriving the genuine low-resolution 
Maxwell theory, it is instructive to again change the no- 
tation to it tot and g tot by adding up t he ma terial and 
field contributions, as we did in section II A 3. The new 



theory remains closed and complete, and is given by the 
Maxwell equations (||^), the continuity equations (p3) 
for mass and entropy, and Eqs (|l]) for it tot and g , 
where the fluxes 



Q 1 



(TS + £,aPa +V-g 



ce x b, 



(33) 



rrtot 



= (-u tm +Ts + i aPa +vg M +e z + h z )5 lk 
+5 4 M ^fe - e 4 e fc - bib k (34) 

are obtained by adding the fluxes of Eq ( |l8| ) and (|2(]) to 
that of (pTj). The total momentum density is still as in 
Eq (p3|) , (p ot = pv + e x b/c, while the energy du tot is 
obtained by adding Eq (|2^) to the differential of the first 
of Eqs (17), yielding 



dtt* 



= Tds + £ a dp a + v ■ dg M + e • de + b ■ db, (35) 
Tds + i a dp a + v ■ dg tot + e • de + b • db. (36) 



Eq (p^) is algebraically identical to Eq (pq), because v ■ 
d(g™~^-g tot )c = -i)'d(exb) = —(vxe)-dh+(vxh)-de, 
where eo = e + (v/c) x b, bo = b — (v/c) x e are the 
respective rest frame fields. For reasons that will become 
clear soon, we shall refer to this set of equations as the 
"2-field" theory. 

In the low-resolution theory of the next chapter, we 
shall be dealing with a dense system containing hidden 
charges and dipole moments. This entails coarse-graining 
operations, of which the most obvious result is the ap- 
pearance of four fields, E,D,B,H, replacing the two 
here. For the remaining of this chapter, however, we 
shall ignore this complication, while considering the less 
obvious, but not less important changes. 
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Of these, the most relevant one is that we no longer 
have an explicit expression for the field energy u F ; more- 
over, it is no longer possible to distinguish u F from u M , 
such that the former only depends on the field variables 
and the latter only on the hydrodynamic variables. (We 
may of course expand u F in the low- field limit, and re- 
tain only the lowest order term, say ^B 2 /p. Although 
this resembles the microscopic expression, i& 2 , the coef- 
ficient p is a function of hydrodynamic variables.) Only 
u tot _ u f _|_ u m j-gtajns a clear-cut significance, because 
it is conserved. This is relevant, as the form of u F was an 
important input of the electromagn etic pa rt of our the- 
ory, cf the last paragraph of section II A 2 . On the other 
hand, this is not an entirely new circumstance, as we did 
also loose the knowledge about the explicit form of u M 
when going over to the dense, second system, cf Eq (28). 

So we shall proceed in two steps. First, we reconsider 
the hydrodynamic theory of a neutral system to under- 
stand its derivation, especially how the ignorance about 
the energy was overcome. Next, this method is general- 
ized to include the e, b-ficlds. The result will be a low- 
resolution Maxwell theory — except that two instead of 
four fields are considered, or the 2-field theory. 

In equilibrium, the entropy s is a function of all the 
conserved quantities in the system, because these are 
the only ones that do not change when the system 
is closed. In our case, the variables are u M ,gf l ,p a . 
Given some (usually satisfied) mathematical properties 



Eq ( pq ) to find that we have already made the right choice 
of variables for the "2-filed" theory. What is more, given 
the cogent, one-to-one relationship between the thermo- 
dynamics and hydrodynamics, the associated set of dif- 
ferential equations also possesses the right structure of 
the low-resolution theory. This is fortunate, as all there 
is left to do is to find out how the two fields turn into four. 
That is, we may now construct the genuine low-resolution 
theory by requiring it to reduce to the present one in the 
2-ficld limit, E, D — > e and H, B — * b, (implying also 
E , D e and H , B -> b ). 

One point remains to be settled: Given the algebraic 
identity, 



= Ts + i a p a + Vig t 



e • e + b • b, 



(39) 



the "2-field" theory must be given only in terms of the 
quantities in this equation. Unfortunately, the two fluxes, 
Eqs d|| fl, do not seem to conform to this rule and 
contain g . This is alright at present, because g M = 
pv is given in terms of two quantities that do appear 
in Eq (|39"|), but it seems to spell disaster for the low- 
resolution theory, because g M is without a clear meaning 
there. 

Similar to the case for it tot , there is no unambiguous 
division of g tot into g M and g F in systems with hid- 
den charges (although it is being attempted over and 
over again with predictably rather contradictory results). 
This is because g M lacks any physical significance: First, 



of .f£, gf, Pa ), we may rewrite this as u-(s, p, g?) t or g u fa devoid of ^ significanc ^ ^ the cc T nserved materia i 



Eq (|29|). [Note that, at this stage, we do not yet know 
the relation between gf 1 and v = du M /dgf 1 .] 

As remarked in the last section, the hydrodynamic 
equations are given in terms of the thermodynamic vari- 
ables and the associated conjugate variables. This is 
closely related to the algebraic identity given by the time 
derivative of Eq (53) , 



u =Ts + t, a p a + Vig, 



M 



(37) 



which the equations of motion 1x24 



must satisfy, 



V • Q M = TV • (sv) + £ Q V • ( Pa v) + WiVj-Ilg. (38) 

This is a very confining expression, and although not 
all that obvious, does uniquely determine the fluxes Q M 
and with the help of some auxiliary thermodynamic 
considerations. (This is frequently referred to as the hy- 
drodynamic standard procedure.) So the hydrodynamics 
is indeed determined once the thermodynamics is; there 
is an one-to-one correspondence between them. 

Given the explicit form of Q M , the equality g M = pv 
may now be established via g M = Q M /c 2 . Had Q M 
been different for some reasons, so would g M , but the 
hydrodynamic theory - in terms of gf 1 and v = du M / dgf 1 
independently - would remain formally unchanged. 

Including charges and fields, the conserved quantities 
are u tot , p a and g tot . So the entropy will depend on 
them and the two fields, s(u tot , p a , g tot , e, b). Again, this 
is equivalent to u tot (s, p a , g tot , e, b) - compare this to 



momentum density of a neutral system. Second, without 
a uniquely defined u M and the associated energy current 
Q M , neither is g M = Q M /c 2 useful as a work-around. 
(Not to mention that this equality is invalid, since £ M 
is no longer conserved). The same holds for g F . Never- 
theless, one clear-cut definition of g M does survive: Be- 
cause of the algebraic equivalence between Eqs ( |35| ) and 
(ptf), we may start from the second expression to obtain 
the first, defining g M = g tot — ex b/c as a shorthand. 
Clearly, this guarantee that the two fluxes, Eqs (33, p4| ), 
are still given by the quantities of Eq ( |39"1 ) , even if g^is 
no longer given by pv. 

Finally, a few remarks about the pressure. While T 



and £, as indicated by Eq (36), remain well defined as 
conjugate variables, and hence retain their place in the 
low-resolution theory of the next section, the pressure 
does not. [This is the reason we have chosen to elim- 
inate it from Eqs ( |32"| , |34| ) .] One might think of tak- 
ing P = -dj ul ot d 3 rfdV, instead of Eq @), only this 
derivative is ambiguous and ill-defined: A dense system 
in the presence of field is anisotropic, so the energy not 
only depends on the volume, as implied by this definition, 
but also on the shape. And the appropriate quantity to 
deal with here is not the pressure, but the stress tensor. 

In isotropic liquids, the pressure encompasses many 
concepts that we find convenient, even intuitive: as the 
surface force density, as the momentum current, as a 
quantity that is continuous across interfaces, and as a 
function of T and p. Hence there is widespread reluc- 
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tance to abandon the pressure at finite fields. Unfortu- 
nately, though there are numerous ways to generalize the 
pressure that will preserve some of these properties, none 
covers all. So one may either define many different pres- 
sures - an approach we eschew as it requires exceeding 
care and has led to considerable confusion in the litera- 
ture - or face up to the Maxwell stress tensor, as we shall 
do here. 



C. Hydrodynamic Maxwell Theory 



Poyntin g ve ctor is obtained by inserting the equations of 
motion (|l], |J, 0); for v = 0, into Eq (||), resulting in 
u ot = - V • Q ot , with Q* ot = cE x H . 

Knowing the rest frame quantities, it is not difficult 
in principle to deduce the corresponding ones for the lab 
frame of a moving liquid, v ^ 0, although this does entail 
some algebra. Fortunately, having established the 2-field 
limit, we have prepared ourselves a shortcut. (The re- 
sults of the dilute systems of course d isplay the correct 
transformation behavior.) First, Eqs ( J36| , f42| ) leave little 
choice other than 



The theory we are going to consider, and the expres- 
sions obtained hereby, are fairly general: They are valid 
for arbitrarily strong fields and nonlinear constitutive re- 
lations; the medium may be moving, and the electromag- 
netic field may depend on time. But we do assume the 
validity of local equilibrium, which needs the character- 
istic time r to be established. Therefore, all frequencies 
are confined to lot <§; 1. (We emphasize that the theory 
is only complete after the dissipative terms have been in- 
cluded. Although the terms given here retain their form 
off global equilibrium, at finite frequencies, we need to 
exercise caution when applying them alone.) 

To derive the general theory, we need, as in the pre- 
vious two cases, the input of the Maxwell equations and 
some information on the energy content of the field. The 
macroscopic Maxwell equations are, 

V • D = p e , V ■ B = 0, (40) 

D = cVxH-p e v, B = -cV x E. (41) 

Taking D, B as the variables, the first two are constraints, 
and the next two their equations of motion. If E and H 
are not specified, Eqs ( [4l|) do not contain any informa- 
tion beyond charge conservation, as they may actually be 
derived from p e + V • (p e v) — and Eqs (^). Especially, 
no restriction whatever is implied for the two "fluxes" E 
and H. 

Taking the temporal derivative of the second constraint 
yields V • B = 0, implying there is a vector field cE 
such that the divergence-free field B may be written as 
B = V x (— cE). Similarly, deriving the first constraint 
yields V- (D+p e v) = 0, so this divergence-free field must 
relate to a field cH, such that D + p e v — V x (cH). 



The temporal Maxwell equations (41 ) are given content 
only by taking E and H as thermodynamic conjugate 
variables, du tot /<9D and <9it tot /(?B, respectively. More 
precisely, we assert that, in the rest frame, the thermo- 
dynamic energy has the form 



dul ot = Tds + (£ Q ) a d Pa + E - dD 



Hn-dBr 



(42) 



It should not come as a surprise that we would need 
this as an input, since the field dependence of u F was 



also an input, cf the last paragraph of section II A 2 . 
(Some textbooks seem able to "derive" the expression 
E • dD + H ■ dB from the Maxwell equations, but this ev- 
idently contradicts our simple consideration above.) The 



dw tot = Tds + £ a dp a + v ■ dg tot 
+E • dD + H • dB. 



(43) 



[The algebraically minded reader may prefer to deduce 
this expression directly. The appendix in section [v| is 
specifically written for him.] Second, since the rest frame 
fields are E = E + (v/c) x B and H = H (v/c) x D, 
we may rewrite Eq (W3) as 



M 



is the generalization of Eq (pq), with g now 



1-EdD + HdB, (44) 
- tot -DxB/c, (45) 

M 



Eq i_ 

given by Eq (|45|). As discussed, this does not imply a 
field contribution of D x B/c, though it does tells us 
that this new g M replaces the old one in the fluxes of 
Eqs ( f46| , f4il| ) below. Note that all conjugate variables are 
now functions of all thermodynamic variables, eg T is 
a function also of D,B, and conversely, Ho depends on 
s,p. 

The rest frame Poynting vector cEo x Ho and Eq ( |3^ ) 
makes it unambiguous that the energy flux, and hence 
the momentum density, now have the form, 



Q 1 



(Ts + £, aPa + v ■ g M )v + c E x H, (46) 



pv + ~E x H/c, 



(47) 



where [Ts + £, a p a + v ■ g ]/ pc 1 is again neglected 
in g tot . (Being a term of zeroth order in the velocity, 
E x H/c may not be neglected with the same argument.) 

To derive the stress tensor, the last unknown expres- 
sion of the low resolution theory, we again take the tem- 
poral derivative of Eq (|43|), u tot = Ts + £, a p a + • • • and 
insert the respective equations of motion, Eqs (^l],^4], [ll]) 
and Eq (Eq), to arrive at 



vg M )v 



■cExH] 



V ■ [{Ts + £, aPa 
= TV ■ {sv) + £ a V ■ { Pa v 
-E • (cV x H p e v) + cH • V x E 



(48) 



This equation is satisfied by the stress tensor 



tot | rri I £ I _M 



= HIT = g™v k - E.Dk - H t B k 



( _ u tot + Ts + ^ pa + v . g M + E . D + H B ),5 lfe (49) 

Although this calculation only yields Vfcll*^, so one can 
always add a term efc mrl V m ^4i„ (with Aj n arbitrary), the 
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requirement that it must reduce to Eq ( |34| ) in the 2-field 
limit eliminates this ambiguity. 

Much information is contained in Eq (^), and one 
could indeed have derived the complete hydrodynamics 
directly from it , without the detour over the dilute sys- 
tem M. The path chosen in this article highlights the 
physics and minimizes the algebra. [To understand how 
stiflingly restricting Eq (48) is, just try out a few terms 



ik -J 

Consider the rotational invariance of the energy it tot 
to see that ±1^ is symmetric: Rotating the system 
by an infinitesimal angle d#, the scalars are invariant, 
dit tot , ds, dp a = 0, while the vectors change according to 
dg M = g u x d dD = D x d dB = B x d Inserting 

these into Eq (pi) yields, 



v x g 



ExD + HxB = 0, 



(50) 



where the left side is equal to 6^11*°*. (For vanishing 
fluid velocity and linear constitutive relations, v = 0, 
eE = D and /iH = B, the symmetry of the stress tensor 
is obvious, otherwise it is not.) 

Aside from the dissipative terms that we have consis- 
tently neglected, the expressions of this section represent 
the complete and closed hydrodynamic theory of dense 
systems that are charged or subject to external fields. 
The independent variables are: s,p a , g tot , D, B, with u tot 
and g tot given by Eqs (|J) and (|^), respectively. The 
attendant equations of motion are the Maxwell equations 
Eqs ©|l|), and 



p a + V • (p a v) = 0, sH 
i tot + V • Q tot = 0, g\° l 



V • (sv) = 0, 

■f v fc n*°.* = o, 



(51) 
(52) 



where Q tot and 11*°,* are given 



by Eqs (H©. 
In comparison to the considerations in available text- 
books, the main advantage of the above derivation is 
the proof of consistency of the equations of motion with 
general principles: Especially thermodynamics, transfor- 
mation behavior, and conservation laws. At the same 
time, the range of validity of these equations are clarified. 
The advantage does not lie so much in the additional, 
frequency and velocity dependent terms, which are fre- 
quently small - although this will change once dissipative 
terms are included. 



D. Boundary Conditions 

To solve the derived set of hydrodynamic equations, 
boundary conditions are needed. In mathematics, boun- 
dary conditions represent information that is extrinsic to 
the differential equations. They are provided to choose 
one special solution from the manifold of all functions 
satisfying the differential equations. In physics, circum- 
stances are frequently different. An example are the 
Maxwell equations, for which the boundary conditions 



are actually derived by integrating the differential equa- 
tions across the interface, yielding 



AD n , AB n , AE t , AH t = 



(53) 



(provided surface charges and currents are absent). The 
subscripts n and t denote the components normal or tan- 
gential to the interface, and AA = — Ar denotes 
the discontinuity of any quantity A, from left to right of 
the interface, toward positive spatial coordinates. Note 
that since the integration has to take place in a frame in 
which the interface is stationary, we have v n — in all 
the boundary conditions derived in this section. 

The results of the integration are better termed "con- 
necting conditions" , since they connect the behavior of 
the variables on both sides of the interface. This is pos- 
sible because the Maxwell equations are valid also on the 
other side of the boundary. In contrast, mathematical 
boundary conditions make no supposition on the validity 
of the differential equations outside the considered region 
- in a sense, they shield the region from outside influence. 
An analogous inte gra tion of the energy and momentum 
conservation, Eqs (|5^) , yields 

AQ*f = 0, AlC = 0, AIC = a(i?r 1 + K, 1 ). (54) 

The first two connecting conditions do not contain much 
independent information: Assuming v — for both sides 
of the interface, and inserting the expressions for Q„, II trl 
from Eqs (^,^), we find them to be satisfied automat- 
ically in equilibrium if the connecting conditions of the 
fields are: Defining m and ti as unit vectors that are 
normal and tangential to the interface, 11*°* = H^tiUk, 
H t = H t U, E t = EiU, B n = B k n k , D n = D k n k , and 
with H t , E t , B n and D n continuous, cf Eqs (|53|), we 
have A<3*°* = and 



AH*°* 



-A(H t B n + E t D n 



0. 



(55) 



Being a generalization of AP — a{R^ 1 + R^ 1 ), the third 
boundary condition for II'°* = n^n^nfe is rather more 
useful and will be employed below repeatedly, (a > 
is the surface tension, and i?i , i?2 the principle radii of 
curvature. The term a(i? 1 " 1 + R^ 1 ) is the equilibrium 
surface current for the momentum density, a singular sink 
or source for the bulk current.) 

This completes the derivation of the hydrodynamic 
theory. Strictly speaking, this is all we need to predict the 
behavior of polarizable systems, as these partial differen- 
tial equations will deliver unique solutions if a sufficient 
number of appropriate initial and boundary conditions 
are provided. This renders our initial question concern- 
ing the force of electromagnetic origin, although as yet 
unanswered, seemingly irrelevant. However, physics does 
not comprise solely of mathematics, and clear thinking 
about the valid expressions for the electromagnetic force 
will give us considerable heuristic power in prediction 
— without having to solve the set of partial differential 
equations each and every time. 
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III. STRESS, FORCE AND PRESSURE 

To find an expression for the force in the presence of 
fields, let us first remind ourselves of its form without 
fields, a quantity that revolves around the pressure: 



gfVvj, is the 



• Its gradient, VP = sVT + p a V£, a 
bulk force density, the quantity that accelerates a 
neutral volume element in Eq (j32|). VP vanishes 
in equilibrium if we neglect gravitation. 

• The pressure P itself is a surface force density, 
which contributes to force equilibrium for instance 
in hydrostatics. P remains finite in equilibrium. 

Since the Maxwell stress IT*^ is the generalization of the 
pressure for polarizable media, we ex pect i ts gradient to 
be a bulk force density, see section 



III A, and a quan- 



tity that relates to t he str ess itself to be a surface force 
density, see section III D| . Presumably, as in the field 



free case, the bulk force density vanishes in equilibrium 
if gravitation is neglected, while the surface force density 
remains finite. 

Generally speaking, the magnetic ponderomotive force 
is up to five orders of magnitude stronger than the elec- 
tric one. This is connected to the fact that their re- 
spective, easily attainable values are similar in SI units: 
E « 10 7 V/m, H « 10 7 A/m (ie E m 30, H « 10 4 , in 
y/j/ m 3 ) . Also, both susceptibilities are similar in mag- 
nitude, and do not usually exceed 10 , hence we have 
PoH 2 w £oE 2 x 10 5 . In addition to the greater ease and 
safety of handling, this frequently makes magnetic fields 
the preferred ones. In almost all the formulas of this and 



the next chapter (except in sections [II A and HIE), the 
electric terms are completely analogous to the magnetic 
ones, and the former may be obtained from the latter 
simply by employing the replacements 



B -> D. H -> E, M 



(56) 



For linear constitutive relations (abbreviated hereafter 
as lcr), they imply p, —> e, x m — * \. To render the 
formulas simple, we shall therefore usually display only 
the magnetic terms. 



A. The Bulk Force Density 



In this section, we shall calculate V^II^' and look for 
all terms in the momentum conservation that alter the 
acceleration p%iV of a volume element, starting from the 
energy density Eq (44), 



Tds + £ a dp a + v-dg M + E-dD + H-dB. (57) 



As discussed in the last chapter, the subscript a enu- 
merates the conserved densities of the system. If there 
are more than one, a summation over a is implied: 
£ a dp a = £dp for a one-component fluid, and £ a dp Q = 
(,dp+ £idpi for a two-component one such as ferrofluids, 



where pi is the density of magnetic particles, p2 that of 
the fluid, and p = p\+ P2 the total density. (A ferrofluid 
is a suspension of magnetic particles, typically of 10 nm 
diameter. Note that p\ is the mass of particles over the 
total volume. Its variation does not arise from compress- 
ing each particle individually, but from increasing the 
number of the particles in a given volume element. This 
appears to be a recurrent misunderstanding in the fer- 
rofluid literature.) It is useful to introduce a number of 
potentials, 

F(T, p a ,Vi, Bi, Di) = u tot -Ts- gfv h (58) 
F(T, p a , Vi,Hi,Ei) = F — H i B l - E.D,, (59) 
G(T,Z a ,v i ,B i ,D i ) = F-Z a p a , (60) 
G(T,£ a ,Vi,Hi,Ei)=F-Z a p a . (61) 

The gradient of G is 

V fe G = -sV k T - p a V k ^ a - g™V k Vi (62) 
—BiVkHi — DiV ' k E.- L . 

Writing the Maxwell stress, Eq ((49|), as 

n*f = -GS ik + gfv k - HiB k - E t D k , (63) 

we find 

V fc n*f = sV t T + p a V£ a + B k {V t H k - V k Hi) (64) 



+D k (ViE k - V k Ei) - Pe Ei + gfViV k + V k {g^v k ) 

which, in addition to the Maxwell equations (fid]), are to 
be inserted into the momentum conservation, 



VfcDS 4 = -pge z 



(65) 



[We shall from now on include gravitation, where g is 
the acceleration of gravity, and the unit vector e z points 
upward.] The result is 

pf [v + (E x H - D x B)/cp] = p e (E + vx B/c) 



(sVT + Pa Vti a + gfVvj) - pge z . (66) 



Compare this to Eq (|32j) and register the great similar- 
ity — but do also remember that the temperature T and 
chemical potentials £ a are now functions of the fields. All 
terms except p4;V may be taken as various force densi- 
ties: The frequency-dependent, so-called Abraham-force, 
— p4f(E x H — D x B)/cp, not usually a large term, the 
"fake" Lorentz force, p e (E + v x B/c), only significant if 
the system is macroscopically charged, p e ^ 0, and the 
bulk force density 



pbulk 



sVT + p Q V£ Q + gfVv.j) - pge z , (67) 



which includes both the gravitational and the electro- 
magnetic force. Note f bulk = — Vfcll^* —pge z in the 
neutral, stationary limit, for v, V x E, V x H = 0, cf 
Eq (|64|). f bulk is the proper macroscopic, coarse-grained 
force, valid as long as the hydrodynamic, macroscopic 
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Maxwell theory is. Together with the Abraham force, 
it accounts for the difference between the true and the 
"fake" Lorentz force discussed in the Introduction. It 
consists only of thermodynamically well defined quanti- 
ties, either variables or conjugate variables of Eq (57). 
This is no longer the case if we follow a widespread con- 
vention to write f bulk as a sum of the zero-field pressure 
gradient and the ponderomotive (or Kelvin) force, as this 
introduces thermodynamically ill- defined quantities and 
unwelcome ambiguities, cf section [II F 3. In equilibrium, 
we have 



p 1 dpx/dt^, which is 



VT,V£i=0, V£ 



(68) 



Inserting these in Eq (p7|), we find, for stationary fluids 



in equilibrium, f 



bulk 



0. 



The considerations of the last two chapters make abun- 
dantly clear that the conserved momentum density is 
the sum of material and field contribution, g tot = pv+ 
Ex H/c, and that the corresponding flux is the Maxwell 
tensor. Nevertheless, in the context of condensed matter, 
because p«>ExH/c, the second term may usually be 
neglected — and since D x B and E x H are of the same 
order of magnitude, so is the Abraham force. [Taking p 
as 1 g/cm 3 , ii as 1 cm/s, H as 10 4 and E as 30, both 
in yff/W, (ie H = 10 7 A/m, E = 10 7 V/m,) we have 
pvc/EH re 3000.] 



B. The Incompressible Limit 

It is noteworthy that the incompressible limit in fer- 
rofluids does not imply p, Vp = 0, and hence V • v = 0. 
This is because incompressibility means the constancy of 
the two local, actual densities, pm of magnetic particles 
and p-p of the fluid matrix (ie pi = (pm) arid p2 = (pf) ; 
with the averaging taken over a volume containing many 
particles) . Yet since the particles are usually denser than 
the fluid, pm ~ 5pF, an increase of particle concentration 
will also increase the total density p = pi + p2. More 
quantitatively, Pi/pm+ P2I Pf — 1 because Pi/pm is the 
fraction of volume occupied by the particles, and pij Pf 
that occupied by the fluid. Taking pm and p-p as constant 
in the incompressible limit, we have 



dp = 7dpi, 7=l-p F /p 



M • 



(69) 



Inserting this expression into £dp+ £idpi, we have £} n dpi 
with £f = £1 + and the modified equilibrium condi- 
tions 



VT = 0, VCi 



(70) 



Ferrofluids may frequently be approximated as being in 
the incompressible limit considered here, because the 
variations of pf and pm, say as a result of field inho- 
mogeneities, are usually negligible. The particle density 
Pi, on the other hand, may vary greatly, and so does p 
via Eq (p9). The parameter controlling this behavior is 



the osmotic compressibility, k c 
larger by around 6 orders of magnitude than the com- 
pressibility of an ordinary liquid, cf section [V A 



C. The Maxwell Stress in Equilibrium 

We consider the simplifications that occur for the 
Maxwell stress tensor if the system is stationary (v = 0) 
and in equilibrium. [Starting from this section, we shall 
no longer display the electric terms explicitly, E, D = 0. 
And the subscript denoting rest frame quantities will 
also be eliminated.! 



1. Total Equilibrium 

With the equilibrium conditions as given in Eqs (|6S 
G is nonuniform only due to the inhomogcncitics in the 
fields, both gravitational and electromagnetic, cf Eq fl62]). 
This we may utilize to obtain a more handy expression 
for the stress. Inserting Eq (|68f ) in (62) and integrating 
the resulting expression from 1 to 2, two arbitrary points 
in the medium, we have 



A(G-pgz + J^BidHi) = 0, 



(71) 



where A denotes the difference between the two points 
of the quantity behind it, eg AG = Gi — G\ . The first in 
the above equation is from AG = f ± VfeGdr^; the mag- 
netic term is from the integral J x (BiV kHijdrk- Writing 

it as L BidHi, we note that it is in fact a difference of 
purely local quantities, with no reference to the path: 
ll B t dH = J 2 BidHi - f* BidHi = A J BidHi, where 
denotes a (possibly virtual) spot of vanishing field, 
Hi — 0, and 1, 2 respectively the field values at the spa- 
tial points 1,2. The integration is unambiguous because 
it is to be carried out in equilibrium, for constant T, £1 
and £ + gz, or for constant T, £™ +jgz in the incompress- 
ible approximation. To emphasize this, the subscript J E 
is added. With Mi = Bi — Hi, we may also write this 
term as \H 2 + J Eq MidHi. The gravitational term comes 

from integrating the chemical potential, — j x pVfe^drfe = 
gtfpdz = gfv pdz -g f pdz = gA(zp), where zp(z) 
= /* pd^ 

Eq (|71J) states the constancy in equilibrium of the 
quantity in the bracket, call it —K, 

-K=G-pgz+±H 2 + J Eq MidHi. (72) 

This enables us to write the stress tensor, Eq ([||), as 

Ilf - {K + \H 2 + J^MidHi - pgz)S ii - HBj, (73) 

the announced handy expression. As we shall see, it 
is rather useful, and even includes the information con- 
tained in the magnetic Bernoulli equation B, see sec- 



tion [II F 2 . More generally, we may calculate II*° for 
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an arbitrary point within the medium — if we know the 
field everywhere and the value of the constant K (usually 
via the boundary conditions to be discussed below). 

If the system under consideration (say the atmosphere) 
is non-magnetic, Mi = 0, B^ = Hi and u tot = u M + 
\H 2 hold. Inserting these into Eq ((72|) and employing 
Eqs (|30|), we have 



K = pgz + P atm , 



(74) 



with K clearly the pressure at z = 0. If the system 
is magnetic, we may still write G(T, £ Q , Hi) = G(0) + 
G om , where the first denotes G for vanishing H, while 
the second the contribution from the field, 

-G cm = J^BidHi = \H 2 + J^MidHi, (75) 

and identify 

K = pgz - (7(0). (76) 

For obvious reasons, G(0) is frequently referred to as the 
"zero-field pressure". This is unfortunate, as G(0) = G 
— G em is a function of T, - these are the natural 
variables of G, and the value of the integral in Eq j75| ) 
explicitly depends on this choice of variables. On the 
other hand, if the zero-field pressure were a physically 
sensible quantity, it should remain the same whichever 
variables it is taken to depend on. It does not, of course, 
as it is equal to G(0). This renders the associated concept 
of pres sure a mbiguous, ill defined, and best avoided, see 
section III F for more details. 

Combining Eq (|6^) and (75) - ignore gravitation and 
the electric field - we find 



Bi&Hi = BiVHi 



(77) 



For a function of one variable, this relation always holds: 
Defining F = J f(g)dg, we have VF = (dF/dg)Vg = 
/Vg. For a function that depends on more variables, 
this relation holds only if the additional variables are 
spatially constant. In the present case, B does depend on 
additional variables, T, £i, £ + gz. Fortunately, all three 
are constant. [B may also, as in next section, depend on 
constant T, p and p\ — which is the reason the subscript 
in J Eq is eliminated from Eq ([77]).] Considering LCR, we 

find V J B l dH l = V(±/iP 2 ), which is equal to BiVHi if 
p (say as a function of T) is constant. 

Turning our attention now to boundary conditions, 
we find that the expression AII*^ [which occurs in the 
boundary condition Eqs (p4j)] may be written as 



AIC = A(K + J Eq M t dH,. 



-Ml 



(78) 



in equilibrium. Two remar ks: (i) The equality A(±# 2 - 
H n B n ) = A(|Af n 2 + \H} - \Bl) was used to derive 



Eq (|T8j) , where A(P 2 ), A(B 2 ) vanish, cf Eq@. (ii) If 
appropriate, substitute J E with J uni , and p with p, cf 
Eq @ below. 



Inserting Eq (|7^) into the third of Eqs (|54l) , assum- 
ing atmosphere on one side, employing Eq (|7J) while ne- 
glecting the gravitational term pgz of the atmosphere, 
the boundary condition reduces to 

K + f^MidHi + \M 2 n - pgz 

- Patm + aiR^ 1 + R^ 1 ), (79) 

a useful formula that we shall frequently refer to below. 

2. Quasi- Equilibrium 

The considerations of the last section is not confined 
to total equilibrium, and may be generalized to include 
quasi-equilibria. Establishing equilibrium with respect 
to the distribution of magnetic (or electric) particles is a 
slow process, because mesoscopic particles diffuse much 
more slowly than atomic ones. Depending on the field 
gradient, the relevant volume and the size of the par- 
ticles, this may take days, even weeks |8|. For a 
rough estimate, we equate the Stokes with the Kelvin 
force to calculate the velocity v with which a magnetic 
particle moves: 6nr]Rv = (AttR 3 /3)- XMoVP 2 /2. Tak- 
ing the particle radius as R = lOnm, the viscosity as 
i] = 10~ 3 kg/ms, the susceptibility as \ the field as 
B = 0.1T, and the field gradient VB as lT/mm, the 
velocity is around 10 -3 mm/s, and the time the particles 
needs to achieve equilibrium is r k lmm/u w 10 3 s. On 
the other hand, particles 10 2 times larger, of the size of 
lp, lead to velocities 10 4 times larger, and a characteris- 
tic time of 10 _1 s. For time scales much smaller than r, 
after a uniform ferrofluid is brought into contact with an 
inhomogeneous magnetic field, the following conditions 
hold instead of Eqs( 



Vpi = 0, VT = 0, pi + pV£ = ~pge z . (80) 

This is because both the diffusion of heat and the estab- 
lishment of mechanical equilibrium are in comparison fast 
processes. The last of Eq @ states f bulk = 0, cf Eq @, 
which is due to the density p quickly turning slightly 
nonuniform to compensate for the field inhomogeneity. 
Inserting Eq (f(| in ©, we find VG = pge z ~ BiVHi, 
an expression that we can integrate directly, from a point 
1 to a point 2, if we approximate the variables T, p, p\ 
as constant, yielding 



A(G-pgz + J un .B l dH l ) = 0, 



(81) 



where the subscript J ■ indicates that the integral is 
evaluated at constant T, p, and p\. Again calling the 
spatially constant quantity in the bracket — K, the stress 
tensor is 

n*f = (K + \H 2 + S^MidHi - pgz)6 ij - H t B, , (82) 

where the subscript of the integral and the lacking bar 
over p are the only differences to Eq (|73|) . Integrating as 
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here the magnetization while holding constant T and p a 
yields the electromagnetic contribution to the free energy, 

-F cm (T, Pa) = f^BidHi = + J^MidHi, (83) 

so K may be written as 

K = pgz-[F(p)-£ aPa ], (84) 

while Eq J74| ) still holds for nonmagnetic media. 

Note that Mj(T, pi,Hi) is the appropriate function for 
evaluating J un .MidHi, while Mi(T,£™, Hi) is the proper 
one for J Eq MidHi - both in the incompressible approxi- 
mation. The first is given by an experiment that quickly 
measures the magnetization while varying the external 
field in a closed system, such that p\ stays constant - 
though contact with a heat bath is necessary to main- 
tain the temperature. The latter, Mj(T, Hi), is to 
be measured in an open system that is connected to a 
particle reservoir, which itself is not subject to a vary- 
ing field, so its chemical potentials £™ remains constant. 
Increasing the field of the system now to measure the 
magnetization, magnetic particles will enter the system 
from the reservoir, resulting in a stronger magnetization 
than in the closed case. 

Either experiment is sufficient to determine both mag- 
netizations, as one can be calculated from the other. This 
is accomplished by employing the thermodynamic rela- 
tion, 
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(85) 



Pi 



where the inverse osmotic pressure is given as k~ s 
x (d^f 1 /dpi)n ■ Eq ( |85| ) may be verified by combining a 
thermodynamic identity with a Maxwell relation, 
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(86) 
(87) 



[We are assuming that M||H, so all M and H here are 
to be understood as the respective magnitude. In addi- 
tion, T is held constant and p is taken to satisfy the in- 
comprcssibility condition Eq ( |69| ) throughout.] For LCR, 
Eq (85|) reduces to 



x ™(T,e?,p) = X m (T,Pi,p)+[pi 1 ^) — 



dpi 

showing that the difference between the two susceptibili- 
ties is of higher order in the field, and may be neglected if 
one strictly adheres to LCR. On the other hand, this cal- 
culation also shows when the second term may no longer 
be neglected: Assuming for an order of magnitude es- 
timate that pi(dx m /dpi) « x m ~ 1, an d tak ing k" 1 
pa 10 3 Pa (see the calculation in section [VA below), 
we find that a field of H « 300 Oe suffices to render 

&/K- 1 « 1. 



D. Equilibrium Surface Force 

Given a solid, polarizable body, or one that is (al- 
though not polarizable) submerged in ferrofluid, we may 
lift the body off the ground electromagnetically, against 
the gravitational force, or balance it against some elastic 
force exerted by a string or spring. The same applies to 
a vessel containing ferrofluid. All these happen even in 
equilibrium, when the bulk force density, Eq (^), van- 
ishes. This force is 



[±M*+f Eq M k dH k }dA 



(89) 
(90) 



= A <j> [J Eq M t dH t + J Eq M n dB n ] dA, 

with the electric counter terms as usual given by Eq (|5£^ 
A now specifically denotes: internal - external, and the 
two integrals are to be taken along the surface of the 
body, one right inside and the other right outside of it. 
As in the last section, f E is to be substituted by / uni , 
and p by p, where appropriate. 

Before deriving Eq (^) , we shall first consider the ram- 
ifications of this formula in two examples: (i) A magne- 
tizable solid body (or a nonmagnetic vessel containing 
ferrofluid) in atmosphere, (ii) a nonmagnetic solid body 
in ferrofluid. 

In the first example, the external magnetization is zero. 
Employing the Gauss law, we change the surface integral 
of Eq @ into a volume integral: With V J Eq M t dH t = 
M t VH t and V/ Bq M„dS n = M n VB n , cf Eq @, the 
electromagnetic force is 



= J (M t VH t + M n VB n )d\. 



(91) 



Consider a plate with the field gradient normal to its 
surface: If the field is predominantly tangential to its 
surface, the electromagnetic force is J M k 'VH k d 3 r; if the 
field is normal to the plate's surface, it is J M k V£?fcd 3 r. 
So interestingly, the magnetic force density interpolates 
between both forms of the Kelvin force as discussed in 



the introduction - and in greater details in section III F 3 
below - though with the difference that these here are 
exact results, as no assumption was made with respect 
to the constitutive relation, especially the density depen- 
dence of the magnetization (or susceptibility). 

If the geometry is more complicated than in a plate, 
say if the magnetizable body is an ellipsoid, it is less clear 
what the normal and tangential component of M, H and 
B are in the bulk, and Eq ( |9l"| ) appears ambiguous. For- 
tunately, this does not matter, because the integral of 
Eq ( |9l"| ) has (by virtue of the Gauss law) a unique value, 
as long as J M t dH t and J M n dB n are given at the sur- 
face, cf Eq (p0|), and the fields M t , H t , M n , and B n being 
continuous in the bulk. [An analogous situation is given 
by the simple integral J fdx, F' — /, and F(l), F(2) 

fixed. Then fdx = F(2) —F(l) is unique irrespective 
of how / varies in the interval between 1 and 2.] 
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Turning our attention now to the second system, a 
non-polarizable body submerged in ferrofluid, we take 
the internal magnetization as zero and find 



-elm 



- f\\K + J Eq M k dH k ]dA (92) 
j> [J Eq M t dH t + / Eq M n dB n ] dA. (93) 



It may also here be useful to employ the Gauss law for a 
conversion of the surface integral into one over the volume 
of the body. However, as the physical fields in the volume 
are not continuous with their respective surface values, 
virtual fields have to be defined which do. They need not 
be healthy fields, may violate the Maxwell equations, or 
have a non-physical susceptibility. 

We now derive Eqs (^) by first considering a solid 
body submerged in fluid and attached to a string. Both 
the solid and fluid may be magnetizable — though differ- 
ently, with the Maxwell stress 11'°* respectively denoted 
as and IlL. One of the boundary conditions at the 
interface is then 



IT 



n 



(94) 



where crf^ = du tot /duij is the elastic stress tensor of the 
solid, given by deriving the thermodynamic energy with 
respect to the elastic strain uy. In equilibrium, VjOy 
= 0. Being a liquid, <rf • = except where the surface 
integral cuts across the string, which also has elasticity 
we need to account for. Same as Eqs (|54|), this boundary 
condition states the continuity of the total momentum 
current cry + Hy . [The surface tension a is irrelevant for 
a solid surface of given shape and therefore eliminated 

Each of the four terms of Eq (|94|) stands for a sur- 
face force density. So the integral over the closed surface 
of the solid body must yield the equation of force equi- 
librium, between the gravitational, electromagnetic and 
elastic force, 



r. 



'last 
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Tf lm = 0. 



(95) 



The first term of Eq (|9J) vanishes, because § cr^dAi 
= J Vja^d 3 r — 0. The last term yields the elastic force 
exerted by the string, which is given by the normal com- 
ponent, Ff ast = — § a f nn dAi, because the string cannot 
sustain a shear stress, a\ n = 0. Clearly, this implies 



Ff m = 4> All -° l dAj . 



(96) 



Because 



' Al^f dAj 



:dA = §{tit k + riin k )AU 



^ An 



^AlTS*) dA 
Eq (|55|), this reduces to 



and because AU^ 



jdA = §8 ik ^ kj 
,dA = §(mAn™ + 
vanishes identically, cf 
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(97) 



Inserting Eq ( |78| ) for A 11^, the Gauss law is em- 
ployed to find AfKdAi = A fS7 l Kd 3 r = 0. As 
a result, only terms that depend explicitly on fields 
are left. From these, the gravitational force is found 
to be .F grav = -fgzApdA = - fV(Apgz)d 3 r = 
— ge z A fpd 3 r, while the electromagnetic force is found 
to be given by Eq (|S9"|). 

Because Eq (|M[) is valid for any interface, and not 
confined to the considered geometry, so are the formu- 
las Eqs (ff), ©, ©, and @. Clearly, there is an 
electromagnetic surface force density of the given form 
whenever the Maxwell stress is discontinuous, ie when 
the magnetization (or polarization) changes abruptly. 

To confirm this, consider a further example: Ferrofluid 
contained in a nonmagnetic vessel, and the whole system 
hung on a string in atmosphere. It is slightly more com- 
plicated geometry, as there are two interfaces to consider: 
ferrofluid-solid, and solid-air. Each has a boundary con- 
dition of the type of Eq |}J). We have cr|^ = cr^ at the 
solid-air interface, because the Maxwell stress is contin- 
uous there, IEy = IEy. An integration yields, as before, 
§a s nn dA l = §al n dAi = ~Tf ast . The boundary at the 
ferrofluid-solid interface is 11? = Ely + Oy , Integrating 
over the interface, we find /(II?. - E^) (L4j +J r f ast = 0, 
or again that Eq J96|), and therefore Eqs (0) and fl89| ) 
to be valid. 

In the literature H, E|, instead of Eq (|S9|), the electro- 
magnetic force is usually given as 



T cXm = jiH.B, - jBkdHkS^dAj. 



(98) 



This expression may also be obtained by integrating over 



] Eq (@). First note §(p* tj +IIy)d^ = / V, (o«. +D- )d 3 
= 0, because the integrand, as the total momentum flux 
density, has vanishing divergence in equilibrium (neglect, 
for simplicity, gravitation). Writing the force equilib- 
rium as /(Elf • + a\j)dAj = 0, again identifying ^ olast = 
— <f> aljdAj, we find J-^ = — §IL$.dAj by comparison 
with Eq (Bq). Now, inserting Eq (O) for 



Hy, and be- 



cause § KdAj = 0, we obtain Eq 

Although algebraically equivalent, Eqs ( |89|) and ( |98| ) 
may behave quite differently when being evaluated for a 
concrete geometry. For instance, considering the force on 
a plate starting from Eq (p8[), it is not easy to find the 
result of Eq (|9l|) . This is because Eq ( |98| ) is a rather non- 
local expression, which does not even explicitly depend 
on the susceptibility of the magnetizable body. 

The above lengthy discussion should not obscure the 
point that these equilibrium surface forces are no more 
than interpretation and visualization of boundary con- 
ditions, especially Eqs ([781 B4|). In fact, when consider- 



ing the experiments in Chapter IV, we shall usually em- 
ploy the boundary conditions, to solidify the algebraic, 
safe, albeit slightly more tedious approach. But we shall 
also frequently point out how the surface forces consid- 
ered here would have yielded the same results. The lat- 
ter serves to demonstrate that these forces represent a 
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heuristic concept of considerable power, and to give fur- 
ther reassurance, if one is needed, of the validity of the 
formulas in this section. 



E. Thermodynamic Derivation of the Stress 

It is worthwhile to rederive our central result, the 
Maxwell tensor, Eq (|63|), via a more direct path, though 
with a rather narrower focus and ensured range of va- 
lidity. This is done in the form of an expanded Lan- 
dau/Lifshitz consideration, separately for the electric and 
magnetic contributions. We start from the energy ex- 
pression Eq (|57|), or from one of its potentials, Eqs (|5§1 , 
|59| , |60| , |6l| ), and shall take pains in eliminating the flaws 
mentioned in the introduction. 

Deforming an isolated system at given entropy, mass 
and electric charge, the energy change is 



SU = S 



dV = - 



y"(An 



Sn)dA k 



(99) 



where the energy density m*°* is given by Eq ([37]) with 
v = 0. Notations: Sri is the infinitesimal (or virtual) 
displacement of the surface, and dA k = n k dA with dA 
the surface element and n k the surface normal pointing 
outwards. The validity of this equation is connected to 
— Ail'j' being the electromagnetic force density needed 
to deform the system, see Eq (|96|). All gravitational 
terms are neglected in this section. 

If systems are considered for which the external stress 
tensor vanishes, because neither material nor field lies 
outside the considered volume, we may substitute n*°* for 
All*]?* in Eq (|99|). Only these systems will be considered 
below. In the absence of fields, II* * — > PSi k with P 
uniform in thermodynamic equilibrium, so SU reduces 
to -P §5r l dA l = -PSV, or Eq ©. Note that if the 
projection of the surface force 11*°* = H^rnrik is positive 
(ie positive pressure if without field), the volume tends 
to expand, if 11*^* < 0, it tends to contract. 

In simple geometries, if it*°* and n*£*<5ri are uniform, 
Eq ( |99| ) reduces to 



SU = 5{u tot V) = u tot SV + V5u tot = -A k U\fSn. 

(100) 

Since u tot is known, we shall evaluate AkH^Sri while 
taking A k and Sr.i to point in all three directions, per- 
pendicular and parallel to the field, and hereby obtain all 
nine components of 11^- ■ 



1. The Electric Part 

We consider a parallel-plate capacitor that is filled with 
a dielectric medium. Denoting the three linear dimen- 
sions of the capacitor as x,y,z, with x <C y,z, the six 
surfaces S^, S^, Sf- (with the outward pointing nor- 



q=Dyz 





FIG. 1: Metal plates at St and with constant charge, 
dielectric medium sandwiched between them. Displacing 
by Sx or Sz respectively decompresses and shears the medium. 



xz, A z — xy, respectively; and the volume is V — xyz. 
Taking the two metal plates as S x , the electric fields 
E, D are along e^, see Fig 1. (We neglect the small stray 
fields at the edges.) The capacitor is placed in a vacuum, 
so there is neither electric field nor material outside. 

We now successively displace the three surfaces S x , 
Sy , in all three directions, Sri = Sx,Sy,Sz (which 
is why the capacitor has to be finite) , while holding con- 
stant the quantities: entropy sV , masses p a V , and elec- 
tric charges q = ±DA X (fr om V D = p c ). Because of the 
simple geometry Eq (100) holds and will be evaluated. 

First, take the surface to be displaced as S% ■ When 
the displacement is Sx, we have 



SV = A x Sx, Ss/s = Sp a /p a = -Sx/x, SD = 0; (101) 

and we have SV , 5s, Sp a , SD = if the displacement is 
Sy or Sz (impl ying a shear motion of S^). Inserting all 
three into Eq (10C), we obtain 



n*°*fe = (Ts + £ aPa - u tot )Sx, 



yrtot 

yx 



n*°* = o. (102) 



If the surface is and the displacement Sz, we have 
SV — A z Sz and Ss/s, Sp a /p a , SD/D = —Sz/z. If the 
displacement is Sx or Sy, we have SV, Ss, Sp a , SD = 0. 
Hence 



n*°*fc = {Ts + i aPa + E X D X ~ u tot )8z, 



(103) 



and n* 



TTtot 

yz 



0. 



As the directions ej, and e 2 



mal ±e x , ±e a , ±e z ) have the areas A x 



yz, 



are equivalent, we know without repeating the calcula- 
tion that a displacement of S+ yields 11*°* = 11*°* and 
n*°* = n*°* = 0. (The term E X D X is a result of the metal 
plates being squeezed, compressing the surface charges, 
Sq/q = SD/D = —Sz/z. The compressibility of the metal 
is taken to be infinite. Otherwise, it would contribute an 
additional term in the stress tensor that we are not in- 
terested in here.) 

These considerations have yielded all nine components 
of 11*°* for a special coordinate system. Because the stress 
tensor of Eq (|63|), for D,E||e,; and B, v = produces 
exactly these components, it is the correct, coordinate- 
independent expression. 
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This conclusion may appear glib, but is in fact quite 
solid. Consider first a vector: If two vectors are shown 
to be equal in a special coordinate system, we know that 
they remain equal in any other system - as long as we are 
sure that they are indeed vectors. The same also holds 
for tensors. (The careful reader may notice an ambigu- 
ity with respect to the off-diagonal part, as both EiDj. 
and EkDi yields the same nine components derived here. 
Fortunately, there is no difference between these two ex- 
pressions, because E x D = for B, v = 0, cf Eq (p0[).) 
This concludes a thermodynamic derivation of the elec- 
tric part of the Maxwell stress tensor. 

Substituting the dielectric medium with vacuum, 11^.°* , 
11'°* reduce to Te 2 /2, respectively, implying a tendency 
to contract along and expand along By and e z — as 
it should in the considered case: The differently charged 
plates want to come closer, while the charge in each plate 
would like to expand. If the medium is one with negligible 
susceptibilities, 11*°*, 11*°* are given as P =p e 2 /2, so the 
same electric force must now contend with the pressure 
in deforming the system. 

Consider the same capacitor, now held at a constant 
voltage 4>. The modified system must lead to the same 
stress tensor, because the stress tensor is a local expres- 
sion and may depend only on the local field. The cal- 
culation is rather similar, though one needs to replace 
m*°* in Eq JlOC| ) with the potential u = u*°* - E D, 
as the system is no longer electrically isolated. And the 
constraint for E as the new variable is Ex — 4> (which re- 
places DA X — q). Connecting the capacitor in addition 
to a heat bath necessitates the employment of the free 
energy, F = u*°* — Ts ED, and changes the constraint 
(from constant sV) to ST = 0. (F is the potential used 
in |5).) For the explicit calculation please cf the mag- 



netic case below, Eqs (104, 105, 106 
replacement B — > D, H — *• E inn 
result is the same, given by Eq (p3| 



, [107], [L08D ) with the 
emcntcd. The final 



2. The Magnetic Part 

To obtain the magnetic part of the stress tensor, con- 
sider a rod along e x , of square cross section, made of 
a magnetizable material and placed in a vacuum. The 
surfaces Sy , Sf are covered with a sheet of metal that 
carries a current J _L e x . With A x <C A y ,A z , the mag- 
netic field will be essentially along and confined to the 
interior of the rod, see Fig 2. So again, there is neither 
field nor material outside. 

If the system is isolated, the metal needs to be su- 
perconducting to sustain the current, and the constraint 
on the variable B during a deformation is constant flux, 
BA X = $. (Compare this to the isolated electric case 
with DA X — q.) If the current J is held constant by a 
battery, the constraint is Hx — J/c (from V x H = j e /c 
and analogous to Ex = <p), and the attendant potential 
is u = w*°* - H B. 

The calculation of the isolated magnetic case repeats 



'Sz 



w\ 



J = Hx 



Sy 



FIG. 2: A magnetic system with constant current J _L x, 
enforced by an external battery, not shown. Again, it is de- 
formed by displacing say S% along x or z. 



the isolated electric one, with all equations - both the dis- 
played Eqs (|l0|Jl0^Jl0|,|l0|) and the others around these 
- remaining valid after the displacement D — > B, E — > H 
has taken place. (When considering compressional dis- 
placements, the elastic energy of the metal sheet is again 
neglected.) 

We shall now consider the rod being deformed at con- 
stant current and temperature, so the energy needed to 
deform the system is 



S(FV) = FSV + V5F = -A k H$5ri, 
SF = -s5T + £ a Sp a — B <5H, 



(104) 
(105) 



where F = u tot — Ts — H- B. We again successively 
displace the three surfaces 5+, S+, 5+, in all three direc- 
tions, Sri — Sx, Sy, Sz, while holding constant the quanti- 
ties: temperature, masses, and the current, ie under the 
conditions, ST = 0, 5(p a V) = and 5{Hx) = 0. 

The first surface to be displaced is S x . When the 
displacement is Sx, we have 



SV - A x Sx, SH/H - 5 Pa / Po 



-Sx/ x; 



(106) 



and we have SV, Sp a , SH = Q if the displacement is Sy 
or Sz (implying a shear motion of S x ). Inserting these 
into Eq (104), we obtain 



n*°*fo 



v , - (Ts + i aPa - u tot )5x, n*°* = n*°* = 0. (107) 

If the surface is and the displacement Sz, we have 
SV = A z Sz, SH = 0, and Sp a /p a = —Sx/x. If the 
displacement is Sx or Sy, we have SV, Sp a , SH = 0. 
Hence 



TtfgSz = (Ts + £ aPa + H X B X - u tot )Sz, 



(108) 



and 11*5°* = n*°* = 0. Since the directions e y and e 2 are 

1*°* = n*°* and n£* = n*°* = o. 



equivalent, we have II* 
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TTtot 



This consideration has yielded all nine components of 
for a special coordinate system. Because the stress 
tensor of Eq (|63|), for B, H||e x and D = produces these 
components, it is the correct, coordinate-independent ex- 
pression. This concludes the thermodynamic derivation 
of the magnetic part of the Maxwell stress tensor. 

Although we have only employed the potentials F and 
F above, it should be clear by now that we could just 
as well have employed G and G from Eqs (^0 61), as- 
suming that the polarizable medium is connected to par- 
ticle reservoirs, which keep the chemical potentials £ Q 
constant. This changes the constraint from Sp a /p a = 
—5x/x, —Sy/y to simply S^ a — 0. The derived expres- 
sion for the stress tensor remains unchanged. 

As a stand-alone, the thermodynamic consideration of 
the last two sections gives us a fairly clear idea on the 
form of the macroscopic Maxwell stress tensor n^j', in 
equilibrium, for v = 0, and with either the electric or the 
magnetic field present. The most important information 
it withholds is that about g\ ot , without which n-£* is 
not unambiguously defined at finite frequencies. Also, 
g\ ot is needed to complete and close the set of differen- 
tial equations given previously, which alone is capable of 
providing a consistent and comprehensive description of 
polarizable systems. In addition, relying solely on this 
consideration, one would perhaps need to be more care- 
ful with the stray fields, especially when thinking about 
the possibilities of terms such as ViEj,ViHj. 

Rosensweig has also considered and derived the mag- 
netic part of the stress tensor in Q . His result is the same 
as the one here, though his algebra is rather more com- 
plicated. More problematically, one of his basic, starting 
assumptions does not hold up: His geometry is a slab 
with A y ,A x <C A z and current-carrying wires along the 
surfaces Sijf, Sf , see his Fig 4.1. The winding of the wires 
is oblique, the currents flow along ±e y in the two larger 
plates S^, but has a component along ±e x in the two 
narrow side walls - take them to be along ±m, a 
vector in the a;z-plane. Rosensweig maintains that the 
resultant field is uniform and perpendicular to the sur- 
face given by the winding, ie by e y and rh. 

Unfortunately, the field is neither uniform nor mainly 
oblique, rendering large parts of the ensuing considera- 
tion invalid. First the qualitative idea: If the two much 
larger plates Si^ were infinite, the field would be strictly 
parallel to . This basic configuration should not change 
much if the plates are made finite, and supplemented with 
the two narrow side walls Sir - irrespective of the cur- 
rents' direction there. This argument is born out by a 
calculation to superpose the fields from various portions 
of the currents. First, divide all currents along ±m into 
two components, along ±e z and ±e x . Next, combine the 
first with the currents along ±e v , such that the four sec- 
tions of the four surfaces form a closed loop at the same 
x-coordinate. The resultant field of all loops is clearly 
the main one, and strictly along e x . The leftover cur- 
rents are those at Sy along ±e x and their total effect is 
a small dipole field. 



F. Ambiguous Notations 

In this section, we shall consider some frequently em- 
ployed expressions and notations that we shall see are 
rather misleading, and therefore perhaps best avoided. 
On the other hand, many of us have so accustomed 
ourselves to this notation that we tend to think along 
the categories it provides. But even then, or especially 
then, one should welcome the opportunity to realize all 
its pitfalls which, as discussed in the introduction, lie 
mainly in the ambiguity of the zero-field pressure, and 
of J M • dH. The expressions one arrives at are for in- 
stance (P+ \H 2 )8ik — HiBk for the stress, or MiVHi for 
the Kelvin force. Scrutinizing their derivation, we shall 
in addition conclude that these formulas are only valid 
for small susceptibility x m <C 1, or M <C H — a range 
of validity dramatically smaller than taken for granted 
usually. So if x m is of order unity, and M ps H , as is 
frequently the case in strongly polarizable systems such 
as ferrofiuids, these formulas are invalid. 



1. Different Zero Field Pressures 

We shall take an approach here that is somewhat 
broader than usual, and simultaneously work with differ- 
ent potentials. This will lead to rather flagrant appearing 
contradictions, the understanding of which should lend us 
a sharpened view of the adopted notation. To simplify 
the formulas, we shall neglect gravitation, and consider 
the incomp ressib le limit of ferrofiuids, p Q ,£ Q — > /?i,£i n , 
cf section [II B| . (Replacing p\ , with p, £ yields the 
results for a one-component, compressible liquid.) 

We start by separating the energy u tot (s, pi, B), the 
free energy F(T, pi,B), and the potentials F(T, pi,H), 
G(T, £i n , B), and G(T, £| n , H), cf Eqs (|| - |6l]), into their 
zero-field and electromagnetic contributions, 



u tot =u(0)+u em (B), F = F(0) + F em (B), (109) 

with analogously defined G cm (B), F cm (H), G om (H). Ad- 
hering to convention, we write the potentials that are 
functions of B as 

«em, F em , G om = JH • dB = \B 2 - JM ■ dB, (110) 

and the potentials with tilde that are functions of H as 

G cm , F cm = -/B ■ dH = -\H 2 - JM ■ dH. (Ill) 

Note that the respective integral is to be taken at con- 
stant s,pi for it cm , at constant T, p\ for F em , F cm , and 
at constant T, ^ n for G cm , G cm - So u cm , F om , G cm are 
in fact not equal, neither are G em , F clrL . The situa- 
tion is similar to evaluating J Tds when considering a 
Carnot process: Since T = T{s,p), we need to specify 
what the density p does when s varies. [We have al- 
ready introduced — G cm = \H 2 + J Eq MidHi and —F cm — 
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\R 2 + J uni MidHi in Eqs (f75|, |8|), given each a specify- 
ing index, and discussed why the difference between these 
expressions is irrelevant as long as LCR holds, cf Eq (|S8|).] 
To evaluate the stress Hik, Eq (|63|), in terms of these 
potentials, we start with its diagonal part — G, given re- 
spectively as 



- G = {s§- s + Pl £- - l) u tot = P(s, Pl ) + \B 
JdB-(l-s 



■\M 2 



' ds 



G=( Pl £--l)F = P(T, Pl ) 



Iff 2 



f2 



, /dB-a-p^M, 



G=( Pl £--l)F = P(T, Pl ) 



(112) 
(113) 



-/dH.(l- Pl £)M 



G = F(T,el n 



iP 



(114) 
/M-dH, (115) 



where depending on the variables, the zero-field pres- 



sures are P(s lP i) = (s4- + P \ 



dpi 



l)u(0), P(T, Pl ) 



(pijfc - 1)P(0), and P(T,Ci n ) = -G(0), cf Eq 



with 



u(0) 



= 0. These four differing G, and with 



them the attendant stress n*£ = —G8ik — P%Bk cer- 
tainly look disturbi ngly contra dict ory. Yet all must be 
equal. Equating Eq(|n|) with (|TT|), we find P(T,£j n ) = 
P{T, P i) — J pi<9M/<9prdH, forcing the conclusions that 
the zero field pressure P is a quantity that depends on 
the variables chosen, and that their difference is of or- 
der magnetic field squared. This means especially that 
only one of the pressures may be field-independent under 
given circumstances. 

The pressure all textbooks (including || ||]) take to be 
field independent is P(T, p), equivalent to P (T, P i) in our 
case. This is the pressure that one would measure in the 
absence of field, for given temperature and density. The 
gradient of the remaining, field dependent terms, Vj [n'°* 
— P(T, P i)dij], cf Eq (114), are then interpreted as the 
electromagnetic force, referred to as ponderomotive or 
Helmholz force. But all this is obviously only valid if the 
temperature and density are indeed kept constant when 
the field is switched on. Under adiabatic circumstances, 
when s is kept constant instead of T, the pressure P(s, pi) 
is the one that is field independent. Consequently, the 
above Helmholz force is no longer the correct expression 
for the electromagnetic force. Rather, one must take the 
field dependent terms from Eq ( |112| ) instead. 

The ambiguity of P arises from the fact that the de- 
pendent variables change when the field is switched on. 
For instance, the chemical potential, £™ = dF(0)/d P i + 
dF cm /d Pl = £i n (0)+A£l n , changes by the amount A£> n = 
dF em /d P = - f dH • ggj-M. Therefore 

p(T,ei n ) = p(T,c(o)) + || r Aei n 

= P(T,p 1 )-/dH.p 1 ^ r M, 



(116) 



explaining the difference between Eqs (114) and (115). 
The ill behavior of P takes one by surprise, as proper 




thermodynamic expressions do not depend on the vari- 
ables chosen. This stems from the fact that P is not a 
thermodynamically defined quantityin the presence of 
fields: It does not appear in Eq ([44]), as all bona fide 
thermodynamic variables do. Since Eq ( ]49| ) or 
sist only of quantities which do appear in Eq 
with or without field - the stress tensor n*£* given there 
holds for any set of variables. These are good reasons to 
make this the expression of choice. 

Frequently, a further approximation is employed for 
the diagonal part of the stress. In dilute, one-component 
systems, the magnetization is usually proportional to 
the density, or M — P (dM/d P ). Similarly, in magnet- 
ically dilute ferrofluids, we may assume that the mag- 
netization is proportional to the particle d ensit y pi, 
M = P i(dM/d P i). Inserting this into Eqs flllj) and 



(114), we find, respectively, 



-G = P+i(P 2 -M 2 ), -G = P+iP 2 



(117) 



The difference between these two expressions are not due 
to a different set of variables, as P is a function of T and 
pi in both cases. So we are dealing with a different pitfall 
here, one that we shall discuss in details in section III F 3 
below. The gist of it is: When we assume M ~ pi, this 
is meant as an approximation, implying the neglect of 
square and higher order terms ~ P \, P \ ■ ■ ■ ■ But consis- 
tency then dictates that we must neglect all higher order 
terms, including M 2 ~ P \. This implies especially that 
the dilute limit is only given if M -C H and x m <C 1 



hold. Therefore, the term lM 2 in Eqs (117) must be 



neglected, and the popular form for the Maxwell stress 
tensor, liif. = (P + jH 2 )8ik — HiBk, is to be taken with 
a large grain of salt, as it is valid only for M H, and 
quite useless if M « H, or \ m ~ !• 



2. Magnetic Bernoulli Equation 

The magnetic Bernoulli equation by Rosensweig is a 
very useful relation. It has been extensively employed 
in his book |2j, and in the literature on ferrofluids. We 
shall include the variation of concentration, which he did 
not consider, and in addition, free this relation from the 



ambiguous notation of the last section III F 1 , in which 
it is given. The point is, the information in the magnetic 
Bernoulli equ ation is contained in Eq (|7l|). Combine it 
with Eq (115) to yield 



A[P(T,ei n )+p.g^]=0, 



(118) 



an expressions of the magnetic Bernoulli equation if the 
system is in equilibrium with respect to pa rticle dis tribu- 
tion. (Note that of all the jMidH t in Eqs (plMT5|), only 
that in Eq (115) has the same variables as J Eq Af,dPj, 
hence these two expressions cancel each other.) Before 
arriving at total equilibrium, and as long as the con- 
centration is uniform, we may start from Eq ( |§l| ) and 
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combine it with Eq (114). The result is 



A[P(T,pi) 



pgz 



px^fMtdHt] 



0. 



(119) 



Substitute p for p\ to arrive at the expression as given 
by Rosensweig. 

The velocity dependent terms in the original magnetic 
Bernoulli equation have not been included here, because 
considerations of mass currents in ferrofluids need to in- 
clude viscosities. Besides, some of the velocity dependent 
terms in the stress tensor, Eq (fi9f), are missing in 0. 



Kelvin and Helmholtz Force 



As in section III F 1 above, we may separate out the 



zero-field pressure P from the hydrodynamic bulk force 
density of Eq @, / bulk = sVT + p Q V£ Q . (We again 
neglect gravitation, and consider a stationary medium, 
v = 0.) The remaining terms are referred to as the pon- 
deromotive force f p , 

/ bulk = VP(T, Pa )-f p , (120) 
f P = (dF eia /dT)VT - Pa V(dF cm /dp a ) (121) 
= V(F om - p a dF em /dp a ) - H.VB,. (122) 

[The relation VF em = (8F eia /dT) ■ VT + {dF em /dp a ) ■ 
V p a + HiVBi was used for the last equal sign.] 

In this section, we shall only employ the two poten- 
tials F(T, p a , Hi) and F(T, p ai B^j. So the zero field 
pressure depends on temperature and densities, and the 
derived ponderomotive force is valid only under isother- 
mal condition of constant densities, cf the discussion of 
section 



IIIF1 



We now proceed to derive the Kelvin force by in- 
corporating some specific simplifications. For LCR, 
-dF em /dT = \H 2 d X m /dT, and similarly dF em /dp a , so 
f p can be cast as 



f p = V{\H 2 p a d X m ldp a ) - \H 2 V X m , 



(123) 



usually referred to as the Helmholtz force in the litera- 
ture. If the system is magnetically dilute, we have x m ~ 9 
in a one-component gas, and x m ~ p\ for dp = jdpi in 
ferrofluids (with no additional dependence on p). Both 
imply Padx/dpa = X- Then the ponderomotive force 
V(H 2 ), or 



reduces to \x 



f p = M^ff,, M = X m H. 



(124) 



Assuming in addition a static field: V x H = 0, the 
force f p may also be written as (MiVj)H, a form one 
often encounters. 

This seems a satisfactory derivation of the Kelvin force, 
as both the linear constitutive relation and the propor- 
tionality to the density are frequently valid approxima- 
tions. Unfortunately, the obviously different expression 



may be obtained by essentially the same derivation. 
We start by defining a slightly different susceptibility: 
M = x m B. Alt houg h different from the more usual con- 
vention of Eq ( 124 ), the new susceptibility is undoubt- 
edly physically equivalent to the old one, and we have 
no a priori reason to prefer either. Both susceptibili- 
ties are related via x m = 1 — \jp, = x m /(l + X m ), or 
dx m = dx m /(p) 2 - So Eq (|123| ) may be rewritten as 

(126) 



f p = V(\B 2 p a dx m /d Pa ) - \B 2 Vr 



This time, assuming ~x m is proportional to one of the 
densities, we obtain Eq (125). 

Which is t he c orrect one, Eq (124) or (125)? Since 



Eq (123) and (126) are algebraically equivalent, the differ- 
ence must stem from the two assumptions, x m ~ 9 an d 
X m = X m /(1 + X m ) ~ 9- Reviewing the above deriva- 
tions, it is obvious that if one of the two assumptions 
were strictly correct, the other would be wrong, and only 
the associated force expression is applicable. Generically, 
however, both x m or % m are power series of p. Assuming 
either susceptibility is linear in p is only an approxima- 
tion. And consistency dictates that all quadratic terms 
are then to be discarded. This implies (i) in the dilute 
approximation, both susceptibilities may only retain the 
term linear in p and are therefore equal; (ii) all other 
terms quadratic in the two susceptibilities are also to be 
discarded, as (x m ) 2 , (x m ) 2 ~ 9 2 —> 0. We conclude: the 
Kelvin force is only valid for X m <g; 1 or M <C H, to lin- 
ear order in x m or M. But then Mi'VBi and M, V-Hj are 
equally valid. The same of course also holds for P^Di 
and Pi'VEi in the electric case. 

The following example aims to illustrate this conclu- 
Consider a thin slab of ferrofluid exposed 
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sion 

to a homogeneous, external magnetic field B, oriented 
normal to the slab. An enforced temperature gradi- 
ent within the slab ensures an inhomogeneous suscep- 
tibility, x m (T)- The internal B- field is uniform, but 
not the internal ff-field, as H = B/{1 + X m )- The 
ponderomotive force of Eq ( |125[ ) is zer o in the given 
circumstance, but not that of Eq (124), which yields 
MiVHi = -x m -ff 2 (Vx m )/(l + X m ), an apparent con- 
tradiction. 

A proper analysis of the force should start from the 
simple expression of Eq j67|). If not, the contribution 
of the zero-field pressure, P(T,pi), in the presence of a 
temperature and concentration gradient, needs to be in- 
cluded. In addition, one must include higher order terms 
in the density, take x = otp + (3p 2 , and x = X./0- ~ x) ~ 
ap+(B+a 2 )p 2 . Inserting these respectively into Eq ( |l23| ) 
and (126), we find, for both cases and a constant B-field, 



f p = \B 2 V[(5p 2 \. 



(127) 



MiVBi 



(125) 



Assuming x ~ 9 or x ~ 9 to hold strictly is equivalent 
to taking (3 — or (3 — —a 2 , resulting respectively in 
f p = and f p = — x-B 2 Vx, or equivalently, Mi'VBi or 

Next we go on to consider nonlinear constitutive rela- 
tions, and shall convince ourselves that the two Kelvin 
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force expressions Eqs (124 125) remain valid - though 
(in a generalization of the above conclusion) only to lin- 
ear orde r in the mag netiz ation and polarization. We in- 
sert Eq ( |110| ) in Eq (122), and equate —dF em /dp a with 
f(dMi/dp a )dBi, to yield 

f p = M t VB t + V / {padMi/dpc - Mi)ABi. (128) 



Assuming that p a dMi/dp a = Mi at constant B - equiv- 
alent to p a dx m /9p a = X m f° r LCR — the integral van- 
ishes, and we retain Eq ( 125 ). If we start from _F cm = 
— / B • dH to evaluate 

f P = (dF em /dT)VT - Pa V(dF em /dp a ) (129) 
= V[F em - Pa dF em /dp a } + BiVHi (130) 
= M t VHi + VJ {p a dM t /d Pa - MjdHt, (131) 

we find f p = M, V-ff, if p a dMi/dp a = Mj for constant 
H. 



IV. EXPERIMENTS 

Having been derived from thermodynamics and con- 
servation of total energy, momentum and angular mo- 
mentum, the theory presented in the preceding chapters 
is fairly general, and valid for all magnetizable and polar- 
izable liquids, from single-component paramagnetic fluid 
to ferrofluids, and for their respective electric counter- 
parts. In the case of ferrofluids, which are suspensions 
of ferromagnetic particles, one is tempted to think that 
the properties of the particles are important, especially 
their orientation (given by the magnetic moment), and 
their internal angular momentum. This is indeed true 
if a high-resolution, mesoscopic account of the system is 
the prescribed goal. However, on a coarser scale, with 
many particles per grain and relevant for most experi- 
ments, the theory derived above is quite adequate, even 
uniquely appropriate for being not unnecessarily detailed 
and complicated. 

As in the last chapter, we shall only display the mag- 
netic terms, accounting for magnetic effects, as the anal- 
ogous electric ones are easily obtained via Eqs (|56|). 



A. Field Induced Variation in Densities 

We consider the change in densities, Ap and Api, of 
a magnetizable liquid in equilibrium, from a region of 
low (or no) field to one of high field. In equilibrium, 
Eqs ( pq) holds. For a one-component liquid, they reduce 
to VT = and 



dp 



Vp 

H 



dH 



VH = -ge z 



(132) 



where the chemical potential has been taken as a func- 
tion of p, T and the field magnitude H . It is useful to 



rewrite the two thermodynamic derivatives, the second as 
[d£/dH] p = -[dB/dp] H = -[dM/dp] H = -H[d X m /dp], 
where the last equal sign is valid only for LCR. The first 
may be approximated: Without field, the inverse isother- 
mal compressibility p 2 [d£ y /dp]T = p[dP/dp]T = n? 1 is 
usually a large enough quantity that one may neglect the 
field related corrections, which is — ^(pH) 2 (d 2 x m /dp 2 ) 
for LCR. Employing this approximation and assuming 
LCR, we integrate Eq (132) to yield the variation in den- 
sity in response to the gradient of magnetic field strength 
and the gravitational potential, 



Ap = p 2 KT A(iif 2 ax m /5p- 



(133) 



where the boldfaced A denotes (as before) the differ- 
ence of any quantity at two different points in the liquid, 
A A = A2 — A\. The electric terms may be obtained, as 
usual, according to Eqs (^6|). (Frequently, this effect - 
referred to as the electro- or magnetostriction — is cal- 



culated using the ponderomotive force of section HI F S . 
The above calculation shows that there is no difficulty at 
all to avoid the ponderomotive force, and the ambiguity 
associated with its notation.) Electrostriction has been 
verified |]l3| , using the refractive index to measure the 
density change. 

If the fluid has two components - such as when it 
is a solution or suspension - we need to consider both 
chemical potentials, V£x = and V£ = — ge z . More 
con venien tly, however, in the incompressible limit of sec- 
tion IIIB, we may consider Eqs ([70|) alone, 



dpi 



Vpi 

H 



dH 



VH 



in 



V Pl = [d^/dp!]^ ({dM/dp^H - 91 e z 



(134) 
(135) 



(Note that dM/dpi is connected to dS^/dH via a 
Maxwell relation.) To our knowledge, field-induced equi- 
librium variation in the solute or particle density has not 
yet been measured in any two-component liquids. This 
is unfortunate especially in ferrofluids, where the vari- 
ation in particle density should be rather pronounced. 
Quantitatively, this experiment yields the thermody- 
namic derivative d^/dpi — [(pi) k os ]~ , with k os the 
osmotic compressibility. Since this is a diagonal deriva- 
tive, its significance in characterizing the ferrofluid ranks 
with that of the compressibility, specific heat and mag- 
netic susceptibility. 

L et us estimate the magnitude of this effect: Writing 
Eq (jT3|) as Api/p! - Kos (p 1 ax m /9pi)A( Mo ^ 2 ), we 
approximate (pidx m /dpi) ~ 1, estimate k os « 10 _3 /Pa, 
and find A Pl / Pl w 0.1 for B = 10~ 2 T. [The value for 
k os is obtained by considering a ferrofluid with 10% of 
its volume occupied by magnetic particles of the radius 
r = lOnm, so the particle density is rti = 0.1/(47rr 3 /3). 
Approximating these particles as ideal gas, the inverse 
osmotic compressibility k^. 1 is equal to the osmotic pres- 
sure, P os = nik^T, so k os = 10" 3 / Pa if T = 300K.] 
Compare this with tiny change of the total density, 
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Ap/p = 5-10 8 at the same field, cf Eq (133), a result 
of the small total compressibility, K = 5 • 10~ iu /Pa. 



B. Current Carrying Vertical Wire 

We consider a vertical wire that goes through a dish 
filled with ferrofluid. Feeding the wire with electric cur- 
rent will drag the ferrofluid toward the wire (located at 
r = and along z in cylindrical coordinates). The sur- 
face of the ferrofluid column is given by z(r), with z at 
which the radius diverges (obtained by extrapolation) as 
the origin, ie z(oo) = 0. The boundary condition Eq ( |79| ) 
is evaluated for two points, z(r) and 0. Because the mag- 
netization vanishes and the curvature radii diverge at 0, 
the attendant result is K = P a tm- Inserting this into the 
boundary condition at z(r), we obtain 



<3> 



a(R\ 1 +R 2 1 ) + pgz, 



(136) 



because M n = for the given geometry. The integral 
is to be evaluated for H = J/(2irr), and at constant 
T, £i if equilibrium has time to be established. For times 
much briefer after the current has been applied, we need 
to substitute p for p, and J uni for J" E , and evaluated 
the integral at constant T, pi, cf the discussion leading 
to Eq (|Bl|). Neglecting the surface tension, a — 0, and 
taking the subscript J . lead to the result give n in M . 

In the spirit of the last paragraph of section [II Ej , we 
remark that Eq ( |136| ) may be considered as an expres- 
sion for force equilibrium, between gravitation, surface 
tension and the magnetic surface force. 

For LCR, Mi = x m Hi, the left hand side reduces to 
^x m H 2 > an d we are left with a quadratic, hyperbolic 
profile of the interface, 8ir 2 p gz = J 2 x m /r 2 if the surface 
tension a is neglected. The effect of a is more impor- 
tant for weak currents, J small. It may be neglected in 
any case for z — > 0, where both curvature radii are large 
enough to be ignored. For z large, one curvature radius is 
simply r, and the other oo. So this part of the ferrofluid 
column is accounted for by 87r 2 (pgz+ a/r) = J 2 x m /r 2 , 
with the term ~ r~ 2 being asymptotically (r — > 0) the 
dominant one. In between, where the actual bend from 
the horizontal to the vertical takes place, both curvature 
radii are finite and need to be included for an under- 
standing of the surface - note, however, that they have 
different signs. 



C. Hydrostatics in the Presence of Fields 

In a system of two connected tubes, with only the sec- 
ond subject to a magnetic field, we expect the ferrofluid 
column to be higher in this tube, as ferrofluid is attracted 
to the region of stronger fields, see Fig 3. To calculate 
the level difference, we employ the boundary condition 
Eq @ for the (flat) liquid -air interface in both tubes. 
Since the field vanishes in the first tube (denoted as 1), 



m 




S-o 



FIG. 3: A U-tube filled with ferrofluid with only the left arm 
subject to a field. The pressure P a tm outside the ferrofluid is 
constant. Depending on the orient ation of th e fie ld, the level 
difference Z2 — Z\ is given by Eqs (138) and (J139|) . 



the boundary condition simply states -P a tm = K — pgz\ . 
Inserting this into the boundary condition for tube 2, we 
obtain 



pg(z 2 - Zl ) = / Eq Midffi + \M 2 n . 



(137) 



Note again the force balance between the electromagnetic 
surface force Eq (p9|), operational at the surface in tube 2, 
and the gravitational force from the disparity in height. 
(As in section III C| , the integral is to be evaluated at con- 
stant T and £i if equilibrium has time to be established. 
For much briefer times after the field has been applied, p 
is to be replaced by p, and the integral by J uni , evaluated 
at constant T, p\.) 

If the field is either predominantly tangential to the 
liquid surface (H = Ht, B = Bt and H n ,B n = 0), or 
predominantly normal to the liquid surface (H = H n , 
B = B n and Ht, Bt = 0), we have respectively, 

gp(z2 -zi) = f^Wm = \H 2 x m (138) 
gp(z 2 - Zl ) = f^MidBi = \H 2 X m ^ (139) 

with the last equal signs in both equations only valid for 

LCR. 

Next we consider the quantity that a pressure gauge 
measures in a ferrofluid exposed to a magnetic field. As 
emphasized, it is not the ordinary pressure, yet as it will 
give some value, the question is what this is. Think of 
the gauge as an enclosed volume of air, at the pressure 
-Patm, see Fig 4. One side of this volume is an elas- 
tic membrane, which is displaced if the external stress 
tensor deviates from the internal one. A finite displace- 
ment d stores up the elastic energy kd 2 /2 per unit area 
of the membrane. (Take the membrane to be stiff, ie 
k large and d small, then we need not worry about 
the pressure change inside.) The elastic energy implies 
a force density kd, rendering the boundary condition 
across the membrane as ALI^ = kd, or via Eq J79|), 
K - P atm + J Eq MidHi + \Ml - pgz = kd. (d is taken 
to be positive when the membrane protrudes into the 
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FIG. 4: A vessel filled with ferrofluid is subject to a field. 
The pressure inside the ferrofluid changes with height and 
the field strength, in response to which the elastic membrane 
of a pressure gauge (right botto m) y ields. The corresponding 
displacement d is given by Eq (140). 



gauge.) We have d = as long as the external stress 
tensor is the same as the internal one, ie in the atmo- 
sphere down to the liquid surface of tube 2, and also just 
below the surface - take this as point 2. Moving down 
the liquid column, to an arbitrary point 3, the membrane 
moves the distance d to maintain force equilibrium. Em- 
ploying the above boundary condition for both points, 
and subtracting one from the other, wc find 



kd = A(f Eq MidH, 



(140) 



with A A = A 3 — A 2 . Note AM„ = Mg-M%, where 
is the magnetization at point 3 normal to the pressure 
gauge membrane, and M% the field at point 2 normal to 
the liquid sur face - both components are not necessarily 
parallel. Eq ( |l40| ) again states a force balance, between 
the elastic, magnetic surface and gravitational force. 

The displacement d is the readout of the pressure 
gauge: If the field is uniform, and if the membrane of the 



gauge is parallel to the liquid surface, Eq (140) reduce 
to pg(z2 — Z3) — kd, the zero-field hydrostatic relation; 
otherwise, field contributions abound - even if the pres- 
sure gauge is simply rotated at point 3 in the presence of 
a uniform field. This amply demonstrates the system's 
anisotropy. A further complication is that all fields are 
the actual ones, distorted by the presence of the pressure 
gauge - though this is an effect that may be minimized. 



D. Magnetic O- Rings and Scrap Separation 



The calculation is already given in section [II D| . Bal- 
ancing the gravitation with magnetic force, we have 

^clm + ^.grav = ^ where ^grav = _ y ^ _ ^ 

J" clm is given by Eqs (§|, ^), and especially by Eqs |§, 
|93| ) if the solid is completely nonmagnetic. (Again, de- 
pending on the time scale, the i ntegral J. may be the 
appropriate one, cf section 1IIC.) 

Because scrap separation is an equilibrium phenome- 
non, we may compare energies instead of computing 
forces. This is a simpler and more qualitative approach 
to understand the behavior of polarizable systems. Con- 
sider for instance the fact that ferrofluid is attracted to 
the region of high magnetic field. Take the field B as 
given, the field energy per unit volume is B 2 /2 in vacu- 
um (or air), and B 2 /(2/j.) in fcrrofmids. Since /1 > 1 
in any paramagnets, the second expression is smaller. 
Given the choice, a volume element of ferrofluid will 
therefore always occupy a region with the highest possi- 
ble field, to reduce the field energy. Conversely, a small, 
non-magnetizable object, submerged in ferrofluid, will on 
the other hand tend to occupy the region of lowest field 
strength. If a difference in height is involved, all these 
of course happen only as long as the gain in field energy 
is larger than the loss in gravitational energy levitating 
the object. (If the field H is given instead of B, we need 
to consider F em = F em — HB. And again, it is larger in 
vacuum then in the ferrofluid: -H 2 /2 > -pH 2 /2.) 
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In this section, we address the physics of some tech- 
nical applications: magnetic O-rings, self-levitation and 
scrap separation. Consider scrap separation first. An in- 
homogeneous magnetic field which becomes weaker with 
increasing height may lift non-polarizable bodies sub- 
merged in ferrofluids off the ground, and hold them at 
specific heights which depend on the shape and density 
of the bodies. 



FIG. 5: Self-levitation of a magnetized body submerged in 
ferrofluid. 

Similarly, a permanently magnetized body submerged 
in ferrofluid tends to collect as much liquid around 
itself as possible, in the space occupied by its stray 
field - even at the price of levitating itself off the bot- 
tom, a phenomenon that is usually referred to as "self- 
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(b) 



FIG. 6: Magnetic pressure seal: (a) principle and (b) enlarged 
displacement of the ferrofluid plug. 



levitation" Q . Frequently, the magnetized body consists 
of a periodic array of north and south poles, with peri- 
odicity A, so the stray field extends one or two A into 
the ferrofluid. Levitated approximately that far from the 
bottom, the body will usually have reached its equilib- 
rium position, as no further gain in field energy may be 
achieved by pushing more ferrofluid below it, and levi- 
tating itself yet higher, Fig 5. 

Generally speaking, although energetic considerations 
are a useful heuristic tool, the actual minimization of to- 
tal energy or free energy is often rather cumbersome, as 
the field does change considerably when a volume ele- 
ment of ferrofluid is displaced. For quantitative calcula- 
tions, it is therefore frequently more convenient to con- 
sider boundary conditions involving the Maxwell stress 
tensor, as we have done throughout this review article. 

This also pertains to magnetic O-rings, which may 
be found in most computer hard disk drives and are 
perhaps the most widely deployed ferrofluid device. In 
these, some ferrofluid is positioned as an O-ring between 
a highly permeable rotating shaft and the pole of a per- 
manent magnet, Fig 6. Serving as a pressure seal, it 
enables the rotary component to work in vacuum. In 
contrast to the above circumstances, we are dealing with 
a metastable state here, as it is always energetically more 
favorable to remove the ferrofluid seal first, have the pres- 
sure equalized, and then return the seal to its original 
place at the poles. So the quantity of interest is the low- 
est energy barrier that must be overcome, which in any 
realistic problems is a quantity notoriously difficult to 
find. The correct force balance, on the other hand, is 
easy to write down: 



AAP 



(^Mn+J Eq M k dH k )dA, 



(141) 



expressing the balance between the difference in the elec- 
tromagnetic surface force [cf Eq (]8{|)] on the two free 
surfaces (area A) of the ferrofluid and the difference in 
external pressure. The equation is obtained by consid- 
ering the boundary condition AIT*°^ = 0, for both free 
surfaces, call them 1 and 2. Use Eq (|78|), forget gravi- 
tation, and remember that P is different for 1 and 2 to 
arrive at Eq (pH|), with AP = P (l) - P Q (2). 



The field is strongest in the middle of the O-ring and 
decays towards both ends. If AP were zero, the force 
F clm would be the same on both surfaces, and the fer- 
rofluid stays in the middle of the O-ring. Increasing 
the pressure on the left (surface 1) pushes the ferrofluid 
towards right, such that surface 1 is in the region of 
higher, and surface 2 is in the region of weaker, fields. 
Equilibrium is achieved when the difference in F° lm bal- 
ances AP. The strongest pressure difference maintain- 
able is when one surface is at the region of highest field, 
and the other is field-free. Assuming for simplicity that 
the magnetic field is predominantly tangential, and that 
LCR holds, we have AP = \x™H 2 . With X m ~ E 
H 2 = hqH 2 , and taking H as of order 10 5 A/m, this 
pressure difference is about 10 5 N/m 2 , approximately the 
atmospheric pressure. 



E. Elliptical Deformation of Droplets 

A droplet of ferrofluid exposed to an external magnetic 
field along e y is distorted, cf eg fli|| . If the field is not too 
strong, it turns from a sphere of radius r to an ellipsoid 
of essentially the same volume, y 2 /a 2 + (x 2 + z 2 )/b 2 = 1, 
with a, 6 the semimajor and semiminor axes, b < a, see 
Fig 7. Given a uniform applied field, B cx , the internal 
field of a magnetizable ellipsoid is also uniform, and along 
the applied field, even for nonlinear constitutive relation, 
cf § 8 of Once these facts are assumed, it is fairly 
simple to calculate the distortion, even nonlinearly, as a 
function of the field, because again, only the boundary 
condition, AU%£ = a(R^ 1 + R^ 1 ) from Eqs (|J), needs 
to be evaluated, at the points (x,y,z) — (0, a, 0) and 
(6,0,0). They are respectively 

AE£? = a(2a/6 2 ), AII^ = a(l/6 + b/a 2 ). (142) 



Employing Eq (^9|) to evaluate the respective left side 
of these two equations, we find An*°* - AIT^ = \M 2 , 
because all fields are constant within the ellipsoid, and 
because M n = M at (0,a,0), M n = at (6,0,0). (The 
extent of the droplet is presumed small enough for us to 
neglect the gravitation.) Taking the difference also of the 
right hand side, we arrive at 



\M 2 = (a/r)(27T 2/3 - V~ 1/6 - V 5/G ), 



(143) 



where the parameter r\ is related to the eccentricity e: 



1 — e 2 = (b/a) 2 . (With r 3 = a6 2 , we have a 



rrj 



-1/3 



rrj '.) This result is the same as that given in [15l 



<1 

6: 

Both sides of this equation may be approximated inde- 
pendently: For small eccentricity, the right side reduces 
to (a/r)2e 2 , while the left side reduces to — 1) 2 H 2 
for LCR. In addition, we have H ex /H = 1 + n(/i — 1)), 
where H ex is the applied field, and n = (1 — 2e 2 /5)/3 for 
small eccentricity, cf §8 in H. 
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Assuming LCR, or u 
have 



\{D 2 /e + pv 2 ) with D = eE , we 



u-v g tot = \{D 2 /e - pv 2 ) - v ■ (E X H)/c 

= \{D 2 /e - pv 2 ) - v ■ (D X H)/ce + G(w 2 /c 2 )(145) 



and deduce 



FIG. 7: A freely suspended ferrofluid droplet elon gates along 
the field, the extent of which is given by Eq (143). 



d(u- v-g tot )/dT> 
= (D + ox H/c)/e 



D /e 



(146) 



APPENDIX 



The validity of Eq (|43|) is shown here directly by trans- 
forming the rest frame expression. More specifically, we 
demonstrate d(u — i>gr tot )/9D = Eq, holding s, p a , v, 
and B constant. We start with 



u (D ^T>, Bq-^B) 



\pv 2 . 



(144) 



This pleasingly simple expression is a result of the ac- 
cidental cancellation of the terms from the Galilean- 
Lorentz transformation with that of the Tailor expansion, 



\pv 2 



u(D,B) = u (T> ,B ) + 2v g^ 
= u (D , B ) + 2v • (E x H)/c + \pv 2 

= u (n,B) + ±gv 2 . 



Higher order terms (such as one ~ -D 4 in the energy u) 
do not invalidate this result. The terms in the magnetic 
field behave analogously. 



In Eq ( |145| ), the explicit form of g in the lab frame 
was employed to deduce the lab frame energy, Eq 



from which then the lab frame energy flux, Eq (46), is 



deduced. This may appear as an inconsistency, but is not, 
because with the rest frame expression for g ot given, we 
already know that the term ~ v is from the rest mass. No 
detailed information about the energy flux is necessary 
here. 
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